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ABSTRACT 


This  thesis  investigates  the  fundamentals  of  multiple-input-multiple-output 
(MIMO)  radio  eommunieation  systems  with  spaee-time  eodes.  A  MIMO  system  was  de¬ 
signed  using  the  Alamouti  spaee-time  eode.  The  modulation  teehnique  was  binary 
phase-shift  keying  (BPSK).  Matlab  with  Simulink  was  used  to  simulate  the  design, 
whieh  was  tested  in  both  an  additive  white  Gaussian  noise  (AWGN)  channel  and  in  a 
multipath  fading  channel  with  AWGN.  Theoretical  performance  was  derived  for  both 
channels  and  compared  to  simulated  results.  The  original  receiver  design  was  changed  to 
incorporate  a  maximal-ratio  combiner  (MRC)  receiving  technique  with  channel  state  in¬ 
formation  (CSI).  The  theoretical  performance  for  this  design  was  determined  and  com¬ 
pared  to  simulated  and  published  results. 
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EXECUTIVE  SUMMARY 


In  modern  radio  communications,  the  demand  for  high  speed,  reliable  communi¬ 
cation  together  with  efficient  use  of  the  spectrum  and  power  are  the  prime  technical  crite¬ 
ria  for  communications  engineers.  However,  these  implementations  must  operate  in  a  se¬ 
vere  environment.  They  must  confront  rich  scattering  effects  and  multipath  propagation 
under  limited  power  constraints.  This  not  only  restricts  the  speed  of  communications  but 
the  reliability  as  well.  Multiple-input-multiple-output  (MIMO)  communication  systems 
with  space-time  coding  provide  increased  reliability  and  increased  capacity  without 
bandwidth  expansion.  MIMO  communications  systems  have  application  to  existing  and 
future  military  radio  communications  systems  where  channel  bandwidth  is  restricted  but 
increased  capacity  is  necessary.  If  the  information  rate  needs  to  double  and  must  fit 
within  the  existing  bandwidth,  then  MIMO  communications  systems  can  provide  the  in¬ 
crease  in  capacity. 

This  thesis  investigated  the  fundamentals  of  MIMO  systems  with  space-time  cod¬ 
ing.  A  simulation  of  a  simple  design  was  created  to  demonstrate  and  analyze  perform¬ 
ance.  The  space-time  code  chosen  was  the  Alamouti  scheme  using  binary  phase-shift 
keying  (BPSK)  with  a  two  transmit  and  two  receive  antenna  MIMO  system.  The  system 
design  was  implemented  in  Matlab  with  Simulink.  The  design  was  tested  in  an  additive 
white  Gaussian  noise  (AWGN)  channel  and  then  a  multipath  fading  channel  with 
AWGN.  The  original  design  was  modified  using  a  maximal-ratio  combiner  (MRC)  tech¬ 
nique  and  tested  in  a  multipath  channel  with  AWGN.  The  simulated  performance  results 
are  compared  with  theoretical  analysis  as  well  as  published  results. 

The  results  showed  that  the  designed  MIMO  system  performed  within  expectation 
of  theoretical  analysis  in  AWGN.  It  outperformed  a  single-channel  BPSK  system,  need¬ 
ing  3.0  dB  less  power  for  equal  performance.  However,  the  designed  system  failed  in  a 
multipath  channel.  The  MRC  design  was  incorporated  in  the  receiver.  As  a  result,  the 
new  design  was  able  to  communicate  successfully  over  the  multipath  channel.  The  per¬ 
formance  of  the  new  design  matched  expected  theoretical  performance  as  well  as  pub¬ 
lished  results.  The  designed  MIMO  system  with  the  MRC  receiver  was  able  to  achieve 
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full  diversity  order  with  the  Alamouti  spaee-time  eode  whieh  was  eonsistent  with  pub¬ 
lished  results.  It  was  able  to  aehieve  a  signifioant  advantage  over  a  single-ehannel  BPSK, 
or  single-input-single-output  (SISO)  system,  in  a  multipath  ehannel  with  AWGN. 


I.  INTRODUCTION 


A.  BACKGROUND 

In  modern  radio  communications  the  demand  for  high  speed,  reliable  communica¬ 
tion  together  with  efficient  use  of  the  spectrum  and  power  are  the  prime  technical  criteria 
for  communications  engineers.  For  example,  the  implementation  of  wireless  Internet 
with  video  conferencing  or  other  media  applications  requires  wideband,  high  data  rate 
and  reliable  communication  systems.  However,  these  implementations  must  operate  in  a 
severe  environment.  They  must  confront  rich  scattering  effects  and  multipath  propaga¬ 
tion  under  limited  power  constraints.  This  not  only  restricts  the  speed  of  communica¬ 
tions,  but  the  reliability  as  well.  Added  to  this  is  the  limited  available  spectrum  that  is  al¬ 
located  for  wireless  systems.  Limited  bandwidth  means  limited  data  rates,  or  speed  of 
communications.  What  if  a  system  can  be  implemented  that  provides  reliability  and  in¬ 
creased  capacity  within  limited  bandwidth  constraints?  This  would  give  the  communica¬ 
tions  engineer  more  avenues  for  trade-offs  in  design.  The  demand  for  increased  capacity 
within  limited  bandwidths  is  certainly  relevant  to  existing  and  future  military  communi¬ 
cations  systems.  A  multiple-input-multiple-output  (MIMO)  communication  system 
with  space-time  coding  holds  such  a  promise  for  engineers.  “Recent  research  in  informa¬ 
tion  theory  has  shown  that  large  gains  in  capacity  of  communication  over  wireless  chan¬ 
nels  are  feasible  in  multiple-input-multiple-output  systems.”[l]  A  review  of  MIMO  re¬ 
lated  literature  finds  the  quote  or  variations  of  it  several  times.  The  primary  focus  of  this 
research  was  to  investigate  the  fundamentals  of  MIMO  systems. 


B,  GOALS  AND  METHODOLOGY 

The  goal  of  this  thesis  was  to  examine  MIMO  systems  with  space-time  coding 
and  their  relevance  to  military  communication  systems.  In  order  to  do  so,  a  simulation  of 
a  simple  system  was  created  to  demonstrate  and  analyze  the  design’s  performance.  The 
space-time  code  chosen  here  is  the  Alamouti  scheme  using  binary  phase-shift  keying 
(BPSK)  with  a  two  transmit  antenna  and  two  receive  antenna  MIMO  system  [I].  The 
modulator  and  simple  receiver  design  were  created  and  implemented  in  Matlab  with 
Simulink.  The  design  was  realized  in  equivalent  baseband  form  to  facilitate  comparison 
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with  published  theory.  The  design  was  tested  first  in  additive  white  Gaussian  noise 
(AWGN).  Theoretieal  performanee  was  determined  and  eompared  with  simulated  results 
of  a  single  ehannel  BPSK  system,  or  single-input-single-output  (SISO)  system.  The  de¬ 
sign  was  then  progressed  to  a  multipath  fading  ehannel  with  AWGN.  The  simulated  per¬ 
formanee  results  were  eompared  with  theoretieal  results  as  well  as  published  results. 

C.  BENEFITS  OF  STUDY 

Spaee-time  eoding  in  MIMO  systems  is  new  teehnology.  There  is  a  great  deal  of 
researeh  effort  in  this  teehnology  [1].  The  demand  for  high  speed  reliable  eommuniea- 
tions  systems  in  a  limited  frequeney  speetrum  for  eommereial  use  is  mimieked  by  mili¬ 
tary  applieations.  However,  some  military  applieations  also  need  eapaeity  and  eoneeal- 
ment.  If  inereased  eapaeity  and  performanee  at  a  eertain  power  level  is  aehieved,  then 
the  power  level  may  be  redueed  by  saerifieing  eapaeity.  Nevertheless,  this  may  be  at  a 
eapaeity  that  is  already  aeeeptable.  The  redueed  power  level  may  lead  into  a  region  of 
low  probability  of  deteetion  eommunieations  and  eoneealment.  Notwithstanding,  MIMO 
and  spaee-time  eoding  holds  promise  for  additional  trade-offs  for  the  eommunieations 
engineer. 

D,  THESIS  ORGANIZATION 

This  thesis  is  organized  in  four  remaining  ehapters.  Chapter  II  introduees  the 
fundamentals  of  the  MIMO  model.  It  diseusses  eapaeity,  the  ehannel  and  the  signifieant 
aspeets  of  MIMO  and  spaee-time  eoding.  Chapter  III  introduees  the  design  of  the 
MIMO  eommunieations  system  using  the  Alamouti  spaee-time  eode  [I].  The  transmit¬ 
ter,  reeeiver  and  logie  of  the  reeeiver’s  ereation  will  be  diseussed.  Examination  of  the 
theoretieal  and  simulated  performanee  of  the  MIMO  system  design  is  detailed  in  Chapter 
IV.  It  progresses  through  performanee  in  AWGN  to  performanee  in  a  multipath  fading 
ehannel  with  AWGN.  Chapter  V  reviews  the  results  and  eoneludes  with  reeommenda- 
tions  for  future  study.  Additionally,  there  are  two  appendiees:  Appendix  A  is  a  list  the 
aeronyms,  and  Appendix  B  eontains  the  bloek  parameters  used  in  the  simulations. 
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II.  MULTIPLE-INPUT-MULTIPLE-OUTPUT 


This  chapter  lays  the  groundwork  for  understanding  and  analyzing  MIMO  sys¬ 
tems  in  a  multipath  fading  environment.  The  MIMO  model  is  introduced  and  the  MIMO 
eapaeity  equation  is  derived.  The  eapacity  of  a  MIMO  system  is  greater  than  that  of  a 
SISO  system  [2].  The  inerease  in  eapacity  of  MIMO  systems  is  the  driving  foree  of 
mueh  of  the  researeh  in  this  area.  Next,  an  overview  of  the  multipath  fading  ehannel  is 
presented.  The  signal  will  be  subjeeted  to  multipath  fading,  so  a  baseband  model  is  de- 
seribed  to  faeilitate  analysis.  Then,  diversity  teehniques  and  the  maximal-ratio  eombiner 
(MRC)  are  deseribed.  A  firm  understanding  is  neeessary  for  analysis  of  the  performanee 
of  the  design.  Finally,  space-time  eodes  are  introdueed.  Spaee-time  coding  is  a  mecha¬ 
nism  to  approach  the  capacity  of  a  MIMO  system  [I].  The  analysis  in  this  thesis  is  per¬ 
formed  in  the  baseband  transmission  model.  Thus,  a  definition  of  the  baseband  transmis¬ 
sion  model  is  introduced  first.  This  ehapter  provides  the  fundamentals  needed  to  analyze 
the  design  of  the  MIMO  system  with  the  Alamouti  seheme. 


A,  BASEBAND  TRANSMISSION  MODEL 

The  baseband  transmission  model  is  the  lowpass  equivalent  of  the  bandpass 
transmission  model  [3].  It  eonsists  of  the  eomplex  lowpass  representation  or  baseband 
equivalents  of  bandpass  signals  and  the  bandpass  ehannel  transfer  funetion  as  illustrated 
in  Figure  1.  In  Figure  I,  the  input  bandpass  signal  is  represented  by  x(t)  and  its  Fou¬ 
rier  transform  A(/) .  The  bandpass  system’s  impulse  response  is  represented  by  h{t^ 
and  its  Fourier  transform  //(/).  Lastly,  the  output  of  the  system  is  represented  by  y  (t) 
and  its  Fourier  transform  T(/)  .  Let  the  baseband  equivalent  of  eaeh  of  the  signals  and 
transfer  function  be  distinguished  by  a  tilde,  i.e.,  x{t) ,  h{t)  and  y{t) .  Similarly,  the 
Fourier  transforms  of  the  baseband  equivalent  system  are  represented  by  A(/),  H{f) 
and  Y (/) . 
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Band— Pass 
System 

y(.0 

X(/) 

r(/) 

x(0 

Complex 

Low— Pass 

y(0 

^(/) 

System 

^(/) 

n/) 

Figure  1.  Bandpass  and  equivalent  lowpass  system  (After  Referenee  [3].). 

The  input  signal  x(t)  is  related  to  its  baseband  equivalent  by  [3] 

x(t)  =  (2.1) 

where  Re  [  ]  represents  the  real  part  of  the  expression  in  the  braekets  and  is  the  ear¬ 
ner  frequeney  of  the  bandpass  system.  The  baseband  equivalent  signal  x(t)  ean  be  rep¬ 
resented  by  its  inphase  and  quadrature  eomponents  (t)  and  Xg  (t)  by  [3] 

x(t)  =  x,(t)  +  7Xg(t),  (2.2) 

where 

x^(t)  =  Re[x(t)],  (2.3) 

and 

Xg(t)  =  Im[i(t)] .  (2.4) 

The  Im[  ]  is  the  imaginary  part  of  the  funetion  inside  the  braekets. 

The  magnitude  of  the  frequeney  response  of  the  bandpass  system  (/)  and  that 
of  the  eomplex  lowpass  equivalent  system  fT  (/)  are  illustrated  in  Figure  2.  The  base¬ 
band  equivalent  transfer  funetion  H{f)  can  be  obtained  by  shifting  the  positive  fire- 
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quency  components  of  the  bandpass  transfer  funetion  //(/)  by  to  be  eentered  about 
zero  [3].  The  transfer  funetion  is  assumed  to  flat  over  the  bandwidth  B  . 


Figure  2.  Bandpass  and  equivalent  lowpass  transfer  funetion  (After  Referenee  [3].). 


From  Figure  1,  the  Fourier  transform  of  the  output  signal  of  the  bandpass  system 
7(/)  and  baseband  equivalent  system  is  given  by  [3] 

Y(f)  =  X(f)H(f).  (2.5) 

The  baseband  equivalent  is  given  by 

Y(f)  =  X(f)H(f).  (2.6) 

It  is  assumed  that  the  bandwidth  of  x(t)  is  small  enough  sueh  that  the  magnitude  and 
phase  of  H{f)  can  be  considered  constant  in  this  bandwidth  of  interest;  speeifically, 

H{f)  =  he^\  (2.7) 


where  h  is  a  positive  real  number  and  cp  is  in  the  range  of  [0,2;r] .  Then  is  given 
by 


h{t)  =  . 


(2.8) 
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where  d{t^  is  the  Dirae  delta  function.  The  baseband  equivalent  output  signal  j)(t)  is 
the  convolution  of  x{t)  and  h{t) ,  given  by 

y(t)  =  |  x{x^h{t 

=  1  x{T^he^‘'’d{t -T)dT  (2.9) 

The  bandpass  output  signal  y  (t)  is  related  to  the  its  complex  equivalent  y  (t)  by  [3] 

=  (2.10) 

Like  the  input  signal  x{t) ,  the  baseband  equivalent  output  signal  y{t)  can  be  decom¬ 
posed  into  inphase  and  quadrature  components  by  [3] 

T(0  =  T/(0  +  2>e(0’  (2.11) 

where 

y,(t)  =  Re[y(t)],  (2.12) 

and 

yg(t)  =  Im[y(t)].  (2.13) 

The  bandpass  system  can  be  obtained  from  its  baseband  equivalent.  The  informa¬ 

tion  of  contained  in  x(t)  and  y(t)  are  preserved  in  x(t)  and  y(t)  [3].  In  the  rest  of  this 
thesis,  the  signals  and  systems  are  represented  by  their  baseband  equivalent. 

B,  THE  MIMO  MODEL 

The  MIMO  model  is  described  in  Figure  3.  The  development  of  the  model  dis¬ 
cussed  in  this  section  is  largely  based  on  Reference  [1].  The  model  consists  of  N,  trans¬ 
mit  and  receive  antennas.  From  each  transmit  antenna,  a  signal  Xj ,  where 
7  =  1, 2, 3, . . . ,  ,  is  transmitted.  The  signal  Xj  represents  the  baseband  equivalent  (t) . 
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(The  function  of  t  has  been  dropped  for  convenience.)  Each  signal  Xj  passes  through 
channels  represented  by  h.j .  It  is  assumed  that  the  bandwidth  of  the  signals  Xj  is 
small  enough  that  Equation  (2.9)  holds  true.  The  symbol  h^j  is  a  complex  number  given 
by 

(2.14) 

and  will  be  referred  to  as  the  channel  coefficients.  At  each  receive  antenna,  white  Gaus¬ 
sian  noise  (AWGN)  n. ,  where  /  =  1, 2, 3, . . . ,  ,  is  added  to  the  signal.  The  symbol  n. 

represents  the  AWGN  baseband  equivalent  n.  (t)  [1].  The  received  signal  is  represented 

as  ,  where  /  =  1, 2, 3, . . . ,  ,  and  is  the  baseband  equivalent  of  the  received  signal  i;  (t) . 


At  any  receive  antenna  i ,  the  received  signal  r.  is  the  sum  of  the  transmitted  sig¬ 
nals  through  the  channel  and  the  AWGN  and  is  given  by  [1] 
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(2.15) 


r,  =  +  /z,2X2  + . . .  +  +  n, 

X, 

r,  ,  ,  1  ^2  , 


Thus,  the  system  as  a  whole  can  be  represented  as  a  system  of  linear  equations  given  by 

[1] 


\\ 

^12 

'Kn, 

" 

^2 

= 

^21 

'^2N, 

X2 

+ 

«2 

^(V,l 

^(V,2  ■ 

^N^N, 

(2.16) 


or,  in  matrix  form,  as 


r  =  Hx  +  n  (2.17) 

where  r  represents  an  x  1  column  matrix,  H  is  an  x  matrix,  x  is  an  x  1  col¬ 
umn  matrix,  and  n  is  an  A^^  xl  column  matrix  [1]. 


Let  each  element  in  x  be  a  zero-mean,  independent  identically  distributed  (IID) 
Gaussian  random  variable  [1].  From  information  theory,  the  optimum  distribution  for 
transmitted  signals  is  Gaussian,  thus  they  are  assumed  to  be  Gaussian  random  variables 
[1].  The  autocorrelation  matrix  R„of  the  transmitted  signal  x  is  defined  by  [1,4] 

R„  =  £[xx‘"],  (2.18) 

where  A*^  is  the  complex  conjugate  transpose  of  the  matrix  A  .  Since  x  is  a  vector  of 
zero-mean  IID  Gaussian  random  variables,  they  are  uncorrelated  [4,5].  The  correlation 
matrix  is  diagonal.  The  covariance  of  x  is  given  by  [1,4] 

C„=R„.  (2.19) 


The  total  average  power  F  of  vector  x  is  given  by  [1] 


N, 

N, 

\,—y 

1  1^ 

V  X- 

=y  E 

j\  ' 

ry 

_  V=1 

,/=i 

(2.20) 
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where  tr(  A)  is  the  trace  of  matrix  A ,  defined  as  the  sum  of  its  diagonal  elements.  It  is 

assumed  that  each  transmit  antenna  transmits  equal  power.  The  power  from  any  one  an¬ 
tenna  is  given  by  [1] 


J  A, 


(2.21) 


Then  is 


Rxx=— (2.22) 

A',  ' 

The  symbol  represents  the  identity  matrix  of  size  N  . 

It  was  assumed  in  the  previous  section  that  the  transmitted  signal’s  bandwidth  Xj 

is  narrow  enough  that  the  frequency  response  of  the  channel  is  flat.  Also  we  assume  that 
there  is  no  attenuation  due  to  path  loss  or  amplification  due  to  antenna  gain  except  as 
measured  in  the  channel  matrix  H .  Furthermore,  we  assume  that  the  channel  representa¬ 
tion  matrix  H  is  fixed  during  a  symbol  period,  i.e.,  each  element  of  H  is  a  fixed  coeffi¬ 
cient.  In  the  discussion  of  a  multipath  fading  channel  to  follow,  the  appropriateness  of 
these  assumptions  will  be  made  clear. 

Let  the  normalization  constraints  be  [1] 

,  ,2 

=N,  for/  =  l,2,...A,.  (2.23) 

7=1 

The  sum  of  the  magnitudes  of  each  channel  coefficient  leading  into  any  one  receive  an¬ 
tenna  is  .  If  the  channel  coefficients  are  modeled  as  random  variables,  then  the  nor¬ 
malization  constraint  will  apply  to  the  expected  value  of  the  left  hand  side  of  Equation 
(2.23)  [1].  This  normalization  is  necessary  for  comparison  to  a  SISO  system. 

The  matrix  H  constitutes  the  channel  state  information  (CSI).  It  is  assumed  that 
there  is  perfect  knowledge  of  the  CSI  at  the  receiver.  While  this  assumption  is  never  true, 
it  can  be  approached.  For  example,  each  coefficient  can  be  estimated  by  sending  training 
signals  to  the  receiver  [1]. 
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The  AWGN  at  each  receive  antenna  is  described  by  a  single  column  matrix  n . 

The  elements  of  n  are  zero-mean,  equal  variance,  IID  complex  Gaussian  random  vari¬ 
ables  [1].  As  with  the  transmitted  signals,  the  autocorrelation  matrix  n  is  given  by  [1,4] 

R..=E[i.i.-].  (2.24) 

Similar  to  the  transmit  matrix  x ,  since  each  element  of  the  noise  is  assumed  zero-mean, 
the  covariance  matrix  is  equal  to  the  autocorrelation  matrix.  Each  element  of  n  is  an  IID 
random  variable  with  equal  variance;  thus,  they  are  independent  and  uncorrelated  [5]. 

The  power  of  each  noise  element  is  equal  to  the  variance.  The  autocorrelation  function  is 

[1,4] 

•  (2.25) 

Let  the  average  power  at  each  receive  antenna  in  the  absence  of  noise  be  repre¬ 
sented  by  which  is  given  by  [1,4] 


Pr=^[rX_ 


=  E 

=  E 


Nt  N, 

^kjXj^hlxl 

J=l  k=l 

N,  N, 

'^^hyx.hlxl 

7=1  k=l 


7=1  k=l 


(2.26) 


where  the  channel  coefficients  have  been  assumed  fixed,  i.e.,  not  random  variables.  The 
(  )  represents  the  complex  conjugate.  Substitution  of  Equations  (2.21)  and  (2.23)  into 
Equation  (2.26)  gives 


v=i 


(2.27) 


10 


Then  the  average  signal-to-noise  ratio  (SNR)  at  eaeh  reeeive  antenna,  represented  by  y , 
is  given  by 

7  =  4  •  (2-28) 

It  was  assumed  that  there  was  no  path  loss  or  gain  through  the  ehannel  with  the  exeeption 
of  that  measured  by  the  ehannel  eoeffieients.  The  total  reeeived  power  /^^,beeomes 


P,.,='LP=N,P  =  t'c{R^). 


The  reeeived  signal  autocorrelation  matrix  is 

R,=E[rr'']. 

Using  Equation  (2.17)  in  Equation  (2.30),  we  obtain 

R^^=E[(Hx  +  n)(Hx  +  n) 

=  E1  (Hx  +  n)((Hx)’'^+n 


(2.29) 


(2.30) 


(2.31) 


Expansion  of  Equation  (2.31)  gives 


R  =  E 


=  E 


Hx(Hx)*%Hxn*^  +n(Hx)*^  +nn*^ 


Hx(Hx)’^  +E[Hxn*^]  +  E  (Hx)’^n  +E[nn’'^ 


(2.32) 


since  the  expectation  of  a  sum  of  random  variables  is  the  sum  of  expectations.  The  trans¬ 
position  of  a  product  of  matrices  has  the  identity  [6] 

(Hxf  =x*"H’".  (2.33) 


Then  Equation  (2.32)  can  be  written  as 

R^^  =  E  [Hxx’^H*^  ]  +  E  [Hxn*^  ]  +  E  [x^^H’^n]  +  E  [nn*^  ]  .  (2.34) 
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If  H  is  deterministic  and  x  and  n  are  independent  of  each  other,  then  they  are  also  un¬ 
correlated;  thus, 


ErHxn’^1  =  HE[xlErn’n 

^  ^  L  j  L  J  .2.35) 

E[x’"H’"n]  =  E[x*"]H*"E[n]. 


Since  both  random  variables  x  and  n  have  zero-means.  Equations  (2.35)  becomes 


E[Hxn’^]  =  0 
E[x"H’"n]  =  0, 

and  Equation  (2.34)  reduces  to 

R^^=HE[xx’"]H’^+E[nn’"]. 

Substitution  of  Equations  (2.18)  and  (2.24)  into  Equation  (2.37)  gives 


(2.36) 


(2.37) 


(2.38) 


Einally,  there  are  two  equations  that  define  the  MIMO  model.  Equations  (2.17) 
and  (2.38); 


r  =  Hx-i-n, 

R„„  =HR„.  H’^+R„„. 


(2.39) 


The  first  is  in  terms  of  random  vectors  and  the  second  is  in  terms  of  the  autocorrelation 
matrices.  With  this  information,  the  capacity  of  a  MIMO  system  can  be  derived. 


C.  CAPACITY  OF  MIMO  SYSTEMS 

The  derivation  of  the  capacity  of  MIMO  systems  presented  here  is  roughly  drawn 
from  that  found  in  Reference  [1].  Eirst,  consider  Equation  (2.17) 

r  =  Hx  +  n.  (2.40) 

Any  x  matrix  H  can  be  decomposed  by  the  singular  value  decomposition  (SVD) 
theorem  into  a  product  of  matrices  [4] 

H  =  UDV’^  (2.41) 
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For  reference,  the  matrix  Q  is  defined  by  [1] 


fHH^  whenA^<A^, 

Q  =  i 

whenA^^  < 


(2.42) 


Note  that  the  singular  values  of  H  are  defined  to  be  the  square  roots  of  the  eigenvalues  of 
Q  [1,4].  Let  the  eigenvalues  of  Q  be  represented  by  Z.  where  /  =  1, 2, 3 . . . .  Now,  Equa¬ 
tion  (2.40)  can  be  rewritten  as 

r  =  UDV’^x  +  n.  (2.43) 

The  matrices  D,  U  and  V  have  the  following  properties  [1,4].  Matrix  D  is  an  x  A^^ 
matrix  of  nonnegative,  real,  singular  values  of  H .  Let  the  singular  values  of  H  be  repre¬ 
sented  by  /?. .  Then  the  matrix  D  is  given  by 


A  0 

0  P2 

D=  0  0 

0  0 

0  0 


when  N^<N^, 


(2.44) 


Pj  0  •  •  •  0  •  •  •  0 

0  Pj  • •  •  0  •  •  •  0 

D=  .  .  .  .  .  .  whenA^^<A^,,  (2.45) 

0  0  •  •  •  •  •  •  0 


Pi  0  •••  0 

0  Pj  •  •  •  0 

D  =  .  .  .  .  when  N,  =  =  N.  (2.46) 

_0  0  •••  p^_ 
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Making  the  assumption  that 


P,>P^...>Pn, 


(2.47) 


it  is  seen  that  some  of  the  higher  numbered  p^  ean  be  zero.  If  the  rank  of  H  is  r ,  then 
there  are  r  nonzero  singularities  [4].  The  matrix  U  is  an  x  unitary  matrix  eom- 
prised  of  the  left  singular  vectors  of  H  [1,4] 


U  = 


u. 


(2.48) 


The  left  singular  vectors  of  U  ,  represented  as  Uj,U2 . .  .u^  ,  are  defined  as  the  eigenvec¬ 
tors  of  HH*^  [1,4].  The  matrix  V  is  an  A^^  x  A^^  unitary  matrix  comprised  of  right  singu¬ 
lar  vectors  of  H  given  by 


V  = 


(2.49) 


The  right  singular  vectors  of  V  ,  represented  as  Vj,  V2 . . .  are  defined  as  the  eigenvec¬ 
tors  of  H*^H  [1,4].  Lastly,  a  unitary  matrix  A  has  the  following  properties  [6]; 


A 


(2.50) 


and 

AA’^=AA^'=I.  (2.51) 

Both  matrices  U  and  V  have  these  properties. 

Let  the  following  transformations  take  place  on  r ,  x  and  n  [1]: 

r'  =  U’^r, 

x'  =  V’^x,  (2.52) 

n'  =  U’^n. 
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Then 


r  =  (u-)"r', 

x  =  (v’^)^'x',  (2.53) 

n  =  (u’")^'n'. 

Substitution  of  Equations  (2.53)  into  Equation  (2.43)  gives 

(U*"  )^‘  r '  =  UDV*"  ( V*"  )^‘  x'  +  (U*"  y  n' .  (2.54) 

Multiplying  both  sides  of  this  equation  by  and  using  the  properties  defined  in  Equa¬ 
tion  (2.50)  and  (2.51),  we  obtain 

U’"  (U’"  y  r'  =  U’^UDV’"  ( V*"  y  x'  +  U’"  (u’"  n' 

Ir'  =  IDIx'  +  In'  (2.55) 

r'  =  Dx'  +  n'. 

Therefore,  Equation  (2.55)  is  an  equivalent  representation  of  Equation  (2.40). 

Expanding  Equation  (2.55)  when  N,<N^,  we  get 
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Similarly,  when  N^<  N^,  Equation  (2.55)  expands  to 


0 

0 


0 

0 

0 


(2.57) 


Equations  (2.56)  and  (2.57)  are  deseribed  graphieally  in  Eigure  4  and  Eigure  5. 


Eigure  4.  The  equivalent  MIMO  model  (After  Referenee  [1].). 
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Y '!' 
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N 

'  Y^^ 

Receiver 

^Y 

0 

Figure  5.  The  equivalent  MIMO  model  <  N,  (After  Referenee  [I].)- 


Note  that  the  0  in  Figures  4  and  5  means  that  no  signal  passes  in  that  ehannel.  For  either 
situation,  or  N^<N,,  the  equivalent  reeeived  signal  at  any  reeeive  antenna, 

where  r  is  the  rank  of  H ,  is 

=  yfZiX'i  + nl  where  z  =  1,2,... r .  (2.58) 

The  rank  of  H  ean  be  no  larger  than  the  minimum  value  between  A,  and 
[1,4,6].  Furthermore,  from  Figures  4  and  5  and  Equation  (2.58),  it  ean  be  seen  that  the 
equivalent  reeeived  signal  at  any  reeeive  antenna  r'  is  a  direet  eouple  between  the 

equivalent  transmit  signal  x'  with  an  applied  gain  faetor  ,  a  singular  value  of  H .  No 

longer  are  the  reeeived  signals  a  linear  sum  of  all  transmit  signals  multiplied  by  a  ehannel 
eoeffieient,  as  was  the  ease  in  the  general  MIMO  model. 

The  equivalent  ehannel  deseribes  r  direet  uneoupled  SISO  ehannels.  The  eapae- 
ity  of  the  system  as  a  whole  is  just  the  sum  of  the  individual  eapaeities  [1].  Shannon’s 
eapaeity  equation  is  [7] 

C  =  IFlog2[l  +  SNR],  (2.59) 
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where  C  is  the  eapaeity  in  information  bits/second,  W  is  the  bandwidth  of  the  ehannel  in 
Hertz  and  SNR  is  the  signal-to-noise  ratio.  Thus,  the  capacity  of  the  MIMO  system  is 
given  by 

C  =  rXlog2[l  +  SNR:],  (2.60) 

/=1 


where  SNR'  is  the  signal-to-noise  ratio  for  the  i-th  equivalent  channel.  To  determine  the 
maximum  spectral  efficiency  CjW  ,  it  is  sufficient  to  find  the  equivalent  received  signal 
and  noise  power.  The  equivalent  received  power  in  the  i-th  channel  is  the  power  of  the 
transmitted  signal  multiplied  by  the  associated  eigenvalue.  In  order  to  obtain  the  power 
of  x' ,  the  autocorrelation  matrix  is  found;  that  is, 

R,y=E[x'x"].  (2.61) 


The  substitution  of  Equation  (2.52)  into  Equation  (2.61)  gives 


R  .  =  E 


(2.62) 


Using  the  transposition  identity  defined  in  Equation  (2.33),  we  see  that  Equation  (2.62) 
reduces  to 


R,,=e[v-xx-v] 
=  V*"e[xx*"]v 

=  v’X.v. 


Substitution  of  Equation  (2.22)  in  Equation  (2.63)  gives 


R  =v*^— I  V 

XX  N, 

=  — V%V. 

N,  ' 


(2.63) 


(2.64) 


The  matrix  V  is  an  x  unitary  matrix;  using  the  identity  in  Equation  (2.51),  we  see 
that  Equation  (2.64)  reduces  to 
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(2.65) 


Rx'x'  ^  • 


Similarly,  to  find  the  power  of  the  equivalent  noise  the  autocorrelation  of  n'  is 
used,  which  is  defined  as 


(2.66) 


Once  again,  substituting  the  transformation  defined  in  Equation  (2.52),  we  get 


R  =E 

n  n 


U*"n(u’"n)’ 

=  e[u*W"u] 

=  u*X„u- 

The  substitution  of  Equation  (2.25)  into  Equation  (2.67)  gives 


(2.67) 


R„v=uXXu. 


(2.68) 


The  matrix  U  is  an  x  unitary  matrix;  using  the  identity  in  Equation  (2.51),  we  see 
that  Equation  (2.68)  reduces  to 


(2.69) 


Therefore,  the  SNR  of  the  i-th  equivalent  channel  is 

IP 


SNR'  = 


(2.70) 


Using  Equation  (2.70)  and  (2.60),  we  find  the  capacity  of  the  MIMO  system  to  be 


C  =  W^log 

i=l 

=^iog2n 


1+ 


l=\ 


1  + 


(2.71) 


where  A.  are  the  eigenvalues  of  Q  and  r  is  the  rank  of  H ,  assumed  to  be  given  by 


r  =  min(A^^,A^,) . 


(2.72) 
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The  function  min(A^^,A^J  represents  the  minimum  between  and  N^. 

Ideally,  Equation  (2.71)  should  be  in  terms  of  the  channel  matrix  H .  Let  the  ei¬ 
genvalue  equation  be  defined  by  [1,6] 

(/lI^-Q)y  =  0  where  y;^0.  (2.73) 

The  matrix  Q  is  defined  in  Equation  (2.42). 

fHH*^  whenA^<A^, 

Q  =  ^  ^  (2.74) 

[h’^H  wheniV,  <iV^. 


Equation  (2.73)  has  a  non-trivial  solution  if  X  is  the  root  of  the  characteristic  equation 
[6].  The  characteristic  equation  is 

;7(;i)  =  det(;iI^-Q)  =  0.  (2.75) 

The  eigenvalues  of  a  square  matrix  Q  are  the  roots  of  the  characteristic  equation 
;  therefore,  the  solution  to  can  be  written  as  [6] 


;?  (  ;l)  =  (;l  -  ;Li  )(;i  -  ;l2  )  •  •  •  (;l  -  =  0 


where  are  the  eigenvalues  of  Q  .  Combining  Equation  (2.75)  and  (2.76),  we 

get 

n(i-l,)  =  det(^I^-Q).  (2.77) 

i=l 

This  is  true  for  all  values  of  X  in  the  real  domain  [1,4]. 

Noting  the  similarities  of  the  left  hand  side  of  Equation  (2.77)  to  (2.71)  and  mak¬ 
ing  the  substitution  defined  by  [1] 


(2.78) 
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in  Equation  (2.77),  we  get 


n(-%-l,)  =  det 


n(i+ 

i=\ 


P 

IP 


i.-Q 


V 


)  =  det 


f  p 
I  + - .Q 


(2.79) 


V 


j 


Therefore,  the  equation  for  the  eapaeity  of  a  MIMO  system  beeomes 

IP 


C  =  Wlog,ll 


i=\ 


1  +  - 


r 

^  P  U 

C  =  lElog2 

det 

I.+ — iQ 

1  )) 

1 

(2.80) 


It  is  appropriate  to  diseuss  three  examples  to  illustrate  the  result.  For  the  first  ex¬ 
ample,  let  the  system  be  a  SISO  system,  i.e.,  a  system  with  a  single  transmit  and  a  single 
reeeive  antenna.  From  Equation  (2.59),  the  normalized  eapaeity,  i.e.,  the  maximum  spee- 
tral  effieieney,  beeomes 


f  =  logJl  +  SNR]. 


(2.81) 


Fet  the  SNR  be  10.0  dB.  Thus  the  normalized  eapaeity  is  3.46  bit/s/Hz. 

For  the  same  example  but  this  time  Equation  (2.80)  is  used.  In  this  ease 
Q  =  1 , =1  and  N,=\.  Using  this  information  in  Equation  (2.80),  we  get 


f 

f  p 

C  =  W\og^ 

det 

1 

1  )) 

=  log2 


1  +  - 


cr 


W 

T  =  log3(l+SNR). 


O  J 


(2.82) 


Thus,  Equation  (2.80)  reduees  to  Equation  (2.81). 

For  the  next  example,  eonsider  a  MIMO  system  with  the  same  number  of  transmit 
and  reeeive  antennas,  that  is,  =  N .  Eet  the  ehannel  matrix  H  be  defined  by  [1] 
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Viv  0  0  •••0 

0  y[N  0  •••0 

0  0  0  ■■■  y[N 


(2.83) 


This  is  equivalent  to  N  independent,  uneoupled  SISO  ehannels.  The  matrix  Q  is  de¬ 
fined  by  H*^H  and  given  by, 


Q  = 


N 

0 


0 

N 


0 

0 


0 

0 

N 


=  Mo¬ 


using  Equation  (2.80),  we  find  the  normalized  capacity  to  be 


w 


Substituting  Equation  (2.84)  into  Equation  (2.85),  we  get 

^  f  -\TT-t 


^  p 

det 

I«+ - tQ 

V 

[  Nal  JJ 

- =  10g2 

W 


=  log2 
=  l0g2 


det 

V  V 


NP 
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f 

P 

det 

2 

Iv 

V 

)) 

y 

'(  P  ^ 

det 

Iv 

V 

A 

)) 

(2.84) 


(2.85) 


(2.86) 


Since  the  determinant  of  a  diagonal  matrix  is  the  product  of  the  diagonals  [8],  Equation 
(2.86)  becomes 


-  =  10g2 

W 


V 


(2.87) 


=  A^log- 


1  +  - 

V 
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The  capacity  of  the  channel  using  a  MIMO  system  as  defined  by  matrix  H  in 
Equation  (2.83)  gives  a  capacity  N  times  that  of  the  SISO  case.  This  makes  sense  since 
the  channel  is  equivalent  to  N  independent  SISO  channels.  The  capacity  increases  line¬ 
arly  with  the  number  of  antennas.  If  =  2  ,  then  the  normalized  capacity  for  SNR  of 
10.0  dB  is  6.92  bits/s/Hz. 


Lastly,  consider  the  2x2  matrix  H  given  by 


(2.88) 


This  represents  a  MIMO  system  with  N^  -  N^=2  .  The  channel  coefficients  are  fixed. 
The  matrix  Q  is  defined  by  H*^H  and  given  by 


H^H 


1 - 

1 

1 

1 _ 

- 1 

1 

1 

_ 1 

- 1 

1 

_ 1 

1 - 

1 

1 _ 

(2.89) 


2  0 
0  2 


Then  Equation  (2.80)  becomes 


r 

^  p 

C  =  IElog2 

det 

1 

1  )) 

W 


=  21og2 


1  + 

V 


(2.90) 


The  capacity  of  the  channel  using  this  matrix  H  gives  the  same  result  as  the  pre¬ 
vious  example  when  N  =  2  .  The  capacity  in  this  case  also  increases  linearly  with  the 
number  of  antennas.  If  N  =  2  then  the  normalized  capacity  for  the  same  SNR,  10.0  dB, 
is  6.92  bits/s/Hz.  If  the  Q  matrix  can  be  reduced  to  a  scalar  multiplied  by  the  identity 
matrix,  then  capacity  is  linearly  proportional  to  the  number  of  transmit  antennas,  maxi¬ 
mizing  the  possible  capacity  of  the  system. 
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D.  THE  MULTITPATH  FADING  CHANNEL 

The  worst  channel  for  any  type  of  wireless  digital  communication  system  is  the 
multipath  fading  channel.  Performance  in  this  type  of  channel  is  much  less  than  that  of 
AWGN  alone  [9,10,1 1].  In  this  section,  an  overview  of  the  multipath  fading  channel  is 
presented.  The  MIMO  system  design  will  be  tested  in  this  channel.  Consequently,  a 
model  of  this  channel  is  derived  in  baseband  equivalent  form  for  use  in  the  simulation. 

Firstly,  radio  waves  may  travel  from  a  transmit  antenna  to  a  receive  antenna  by 
several  different  paths.  This  is  called  multipath.  Multipath  propagation  of  the  signal  is 
caused  by  reflectors  and  scatters  present  in  the  physical  channel  [9,10,1 1].  The  signal, 
traveling  via  several  different  paths,  will  be  sensed  at  the  receive  antenna  several  differ¬ 
ent  times  with  different  amplitudes,  phases  and  arrival  times  [9,1 1].  This  results  in  a  re¬ 
ceived  signal  that  can  vary  dramatically  in  amplitude  and  phase.  In  the  frequency  do¬ 
main,  the  spectral  components  of  the  signal  may  be  affected  differently  by  the  channel. 
Thus,  the  frequency  response  of  the  channel  may  not  be  flat  over  the  bandwidth  of  the 
signal.  The  overall  result  of  multipath  propagation  is  time-spreading  of  the  signal  that 
leads  to  intersymbol  interference  [9,10]. 

The  configuration  of  the  reflectors  and  scatters  in  the  physical  channel  may 
change  over  time  as  well.  This  can  be  caused  by  the  reflectors  and  scatters  moving,  or 
more  likely,  the  transmit  or  receive  antennas  moving.  As  a  result,  the  multipath  channel 
is  time-varying.  The  motion  between  the  transmitter  and  receive  antennas  will  cause  a 
Doppler  shift  in  frequency  proportional  to  the  relative  velocity  between  them  [9,10,1 1]. 

The  time-spreading  and  time-varying  aspects  of  the  multipath  fading  channel  can 
be  measured  and  statistically  estimated  [9,10,1 1].  A  maximum  for  the  Doppler  shift  can 
be  determined  as  well.  Knowledge  of  this  information  aids  the  engineer  in  making  key 
design  decisions  for  a  communication  system. 

The  time-spreading  nature  of  the  channel  is  characterize  by  the  coherence  band¬ 
width.  “The  coherence  bandwidth  is  a  statistical  measure  over  which  the  frequency  re¬ 
sponse  of  the  channel  is  considered  fiat”  [11].  The  spectral  components  of  a  signal  that 
fit  in  the  coherence  bandwidth  will  be  passed  with  generally  equal  gain  and  linear  phase. 
If  a  signal  spectrum  fits  within  the  coherence  bandwidth,  then  it  undergoes  fiat  fading. 
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Thus,  the  channel  is  flat  fading  or  frequency-nonselective.  If  the  spectrum  does  not  fit 
within  the  coherence  bandwidth,  then  it  suffers  frequency  selective  fading.  Thus,  the 
channel  is  frequency-selective  [11]. 

The  time-varying  nature  of  the  channel  is  described  by  the  Doppler  spread  and 
coherence  time  [9,10,1 1].  The  coherence  time  is  inversely  proportional  to  the  maximum 
doppler  shift  [9,10,1 1].  The  coherence  time  is  the  statistical  measure  of  the  duration  dur¬ 
ing  which  the  channel  is  essentially  time-invariant.  Thus,  the  coherence  time  gives  an 
indication  of  how  long  the  channel  remains  constant.  If  the  symbol  duration  of  a  signal  is 
less  than  the  coherence  time,  then  it  undergoes  slow  fading.  Thus,  the  channel  is  called 
slowly  fading.  If  the  symbol  duration  is  larger  than  the  coherence  time,  then  the  signal 
suffers  fast  fading  and  the  channel  is  called  fast  fading  [9,10,1 1]. 

The  simulations  were  run  for  a  slow  fading,  frequency-nonselective  channel.  To 
simulate  a  slow  fading  channel,  the  model  of  the  channel  coefficients  is  not  permitted  to 
change  within  a  symbol  period.  This  corresponds  to  the  assumptions  made  concerning 
the  channel  coefficients  in  Sections  A  and  B.  In  order  to  simulate  the  frequency- 
nonselective  or  flat  fading  nature  of  the  channel,  the  channel  coefficients  must  be  distrib¬ 
uted  according  to  a  probability  distribution  function.  Flat  fading  is  generally  simulated 
using  a  Rayleigh  distribution  for  the  magnitude  of  the  channel  coefficients  [9,10,1 1]. 

This  is  consistent  with  a  channel  where  the  fading  is  caused  by  randomly  moving  scatter- 
ers  and  reflectors  or  moving  transmitter  and  receiver.  In  order  to  achieve  the  baseband 
simulation  of  a  Rayleigh  distribution,  two  zero-mean,  independent  Gaussian  random  vari¬ 
ables  are  summed  [11]. 

Let  Z  be  the  complex  sum  of  two  IID,  zero-mean  Gaussian  random  variables 
given  by 

Z  =  X  +  jY  =  He^\  (2.91) 

The  mean  and  variance  of  X  and  Y  are 

(2.92) 

(2.93) 


X  =  T  =  0, 


^2  ^2  1 
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Let  W  be 


r=z^ 


(2.94) 


or,  equivalently 


W  =  X^  +  Y\ 


(2.95) 


The  random  variable  W  has  a  chi-squared  distribution  given  by  [10] 


fw{w)  =  - 


1 


UJ 


w 


,{n/2)-l  -w/2a^ 


i{w), 


(2.96) 


where  the  unit  step  function  is  denoted  by 


!(w)  =  |° 


w  <  0 
w  >  0. 


(2.97) 


The  variable  n  is  defined  as  the  number  of  independent  random  variables  in  the  sum.  In 
this  case  n  =  2  .  The  distribution  then  becomes 


fw{w)=- 


1 


2  '^2/2t-' 

a  2  1 


w 


,{2/2)-l  -wI2ct‘ 


u(w) 


(2.98) 


This  reduces  to 


fwM 


o-"2r(i) 


0  —w/2o  /  \ 

w e  u\wj . 


Since 


r(i)  =  (i-i)!  =  i, 

the  probability  distribution  function  of  W  is  given  by 


(2.99) 


(2.100) 


(2.101) 
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The  magnitude  //  of  Z  in  Equation  (2.91)  is  of  interest  and  is  given  by 

H  =  \Z\  =  ^.  (2.102) 

The  probability  density  funetion  for  W  ean  be  transformed  to  yield  the  probability  den¬ 
sity  funetion  of  H  aeeording  to  [12] 

(2.103) 

where 


/hW  = 


1 


dh  /  dw\ 


fwi^) 


dh  _  1 


(2.104) 


Using  Equation  (2.101)  and  the  substitution  defined  in  Equation  (2.104)  in  Equation 
(2.103),  we  get 

,4  (*)  =  2(A=)'"^e-‘' 

(2.105) 

(J 

The  mean  and  varianee  of  H  are  given  by  [11] 

(2.106) 

Thus,  H  follows  a  Rayleigh  distribution  as  required,  and  the  sum  of  two  zero-mean,  IID 
Gaussian  noise  sourees  will  produee  the  desired  ehannel  model. 
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The  phase  0  of  Z  is  given  by 


Let  A  be 


' 

1  1 

tan 

uJ 

tan  ' 

1  I 

tan  ' 

uJ 

-n 


+  ;r 


for  -oo<x<0,-oo<j;<0 
for  0<x<oo,-oo<y<oo 
for  -  00  <  X  <  0, 0  <  y  <  00. 


A  = 


Y 

X  ■ 


(2.107) 


(2.108) 


For  the  first  ease,  when  -oo  <x<0,-Go<y<0,  the  cumulative  distribution  function  is 
given  by  [5] 

^A{^)  =  f_^£/xY{xx)dydx.  (2.109) 


The  probability  distribution  function  is  given  by  taking  the  derivative  of  F^{a) 


dF^[a)  _  d 
da  da 


By  the  Leibniz  rule,  the  probability  density  function  is  given  by  [5] 

7.4  (« )  =  L  ax)  dx . 


(2.110) 


(2.111) 


The  random  variables  X  and  Y  are  independent;  therefore  their  joint  probability  density 
function  is  given  by 


fxY{^^y)  =  fxi^)fY{y)-  (2.112) 

The  substitutions  of  Equation  (2.1 12)  as  well  as  the  Gaussian  probability  density  function 
into  (2.111)  gives 


/a  («)  = 


1 

4;rcr^ 


{2x)dx 


(2.113) 
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Evaluating  Equation  (2.1 13),  we  obtain 


Ia  («)  = 


2;z'(l  +  a^^ 


0  <  a  <  00  , 


The  random  variable  0  is  given  by 


0  =  tan  '  . 


where 


da 


1  +  fi! 


2  ’ 


and 


=  tan(^  +  ;z-) , 


Using  the  transform  defined  in  (2.103),  we  get 


U(0)  =  \+. 


27i:{\  +  a^^ 


=tan(^) 


Equation  (2.1 18)  simplifies  to 


for-;r<^<-|. 


(2.114) 


(2.115) 


(2.116) 


(2.117) 


(2.118) 


(2.119) 


Eor  the  second  case,  when  0<x<oo,-oo<y<Go,  the  cumulative  distribution  function  is 
given  by  [5] 

Fa  («)  =  £  \_^fxY  {X’  y)  dydx .  (2. 120) 

Eike  the  first  case,  the  probability  distribution  function  is  gained  by  taking  the  derivative 
of  [a] 


dFA{a) 

da 


A. 

da 


f^[x,y)dydx 


(2.121) 
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By  the  Leibniz  rule,  the  probability  density  funetion  is  [5] 


/^(«)  =  Jo  xf^{x,ax)dx. 


(2.122) 


Similar  to  the  first  ease  first  ease,  f^{a)  is 


(2.123) 


In  this  ease  0  is  given  by 


0  =  tan  '  , 


(2.124) 


where 


^  1 
da  \  +  a^ 


(2.125) 


=  tan (6*) . 


(2.126) 


After  transformation  /^{d)  beeomes 


(2.127) 


For  the  last  ease,  when  -co  <x<0,-Go<y<0,  the  eumulative  distribution  funetion  is 


given  by  [5] 


(^)  "  \-Jax^xY  {x,y)dydx . 


(2.128) 


The  probability  distribution  funetion  is  given  by  taking  the  derivative  of  F^{^a) 


dFAa)  <i  r f“  f*  /•  /  \  " 

J  fxY{^^y)dydx  . 


(2.129) 


Onee  again,  by  Leibniz’s  rule,  the  probability  density  funetion  is  [5] 


7.4  («)  =  Jo  -^fxY  ax) dx , 


(2.130) 
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and  f^{a)  becomes 


Ia  («)  = 


27t[\  +  a^^ 


-00  <  a  <  0  , 


In  this  case  0  is  given  by 


where 


0  =  tan  '  (^A)  +  7r , 


da 


1  +  (2 


2  ’ 


and 


=  tan  (6*) 


Finally,  after  transformation  /©(6’)  becomes 


(2.131) 


(2.132) 


(2.133) 


(2.134) 


U(0)  =  ^  for^<0<;r.  (2.135) 

In  1 

Therefore,  the  complex  sum  of  two  zero-mean,  IID  Gaussian  random  variables 
gives  a  magnitude  that  is  Rayleigh  distributed  with  a  phase  that  is  uniformly  distributed 
in  \-n,n\ .  This  models  the  baseband  equivalent  slow  fading,  frequency-nonselective 
channel.  It  will  be  used  to  simulate  a  multipath  fading  channel. 


E,  DIVERSTIY  AND  ERROR  CONTROL  CODING 

As  stated,  the  multipath  fading  channel  is  a  severe  environment  in  which  to  com¬ 
municate.  The  error  rate  performance  suffers  greatly  compared  to  an  AWGN  channel 
[9,10,1 1].  However,  two  techniques  available  to  mitigate  the  effects  of  multipath  fading 
on  performance  are  diversity  and  error  control  coding  [11].  The  possibility  exists  that  a 
MIMO  system  can  use  both  of  these  techniques  in  order  to  improve  performance  [1,13]. 
In  this  section,  diversity  techniques  will  be  investigated.  The  Alamouti  scheme,  the  focus 
of  the  simulations,  incorporates  diversity  techniques  in  its  space-time  code,  so  the  per¬ 
formance  of  diversity  techniques  is  of  great  interest  [1,13]. 

31 


The  concept  of  diversity  is  to  either  transmit  and/or  receive  the  same  information 
more  than  once.  If  one  transmission  suffers  from  a  deep  fade,  then  it  is  possible  that 
other  independent  faded  transmissions  will  not.  The  likelihood  that  all  transmissions  will 
be  faded  below  possible  recovery  is  reduced.  The  receiver,  then,  can  choose  from  the 
best  receptions  [9,10].  There  are  three  methods  of  diversity:  time,  frequency  and  space. 
Time  diversity  transmits  the  same  information  several  times  [9,10].  Frequency  diversity 
transmits  the  same  information  on  several  carriers  [9,10].  Lastly,  space  diversity  em¬ 
ploys  several  antennas  at  the  receiver  in  order  to  receive  the  same  information  several 
times  [9,10].  In  a  MIMO  system,  multiple  antennas  are  employed  so  space  diversity  is 
possible  [1,13].  Also,  frequency  or  time  diversity  is  possible  [1,13].  It  will  be  seen  that 
the  Alamouti  scheme  uses  both  time  and  space  diversity  [1,13].  Therefore,  the  bit  error 
performance  of  time  and  space  diversity  in  a  multipath  fading  channel  will  be  considered 
for  BPSK. 

1,  Maximal-Ratio  Combining 

The  receiver  that  achieves  the  best  performance  in  a  multipath  fading  channel  is 
one  that  uses  maximal-ratio  combining  (MRC).  The  MRC  receiver  uses  channel  estima¬ 
tion,  i.e.,  CSI,  of  the  magnitude  and  phase  of  each  diversity  reception.  The  MRC  receiver 
multiplies  the  received  diversity  reception  with  the  complex  conjugate  of  the  estimation. 
As  a  consequence,  the  phase  is  corrected  for  coherent  detection  and  the  magnitude  of  the 
signal  is  weighted  by  the  signal  strength.  Strongly  received  signals  will  be  weighted 
more  than  weaker  signals  (deep  faded  signals)  in  the  decision  process  [9,10]. 

The  time  diversity  reception  of  one  of  /  =  1, 2, 3, . . .  L  diversity  signals  for  a  BPSK 
system  is  illustrated  in  Figure  6.  Each  received  diversity  signal  r,  (t)  =  y,  (t)  +  (t)  has 

duration  T  seconds.  The  correlation  receiver  equivalent  of  the  matched  fdter  is  depicted. 
Assume  that  the  signal  after  the  correlation  receiver  is  sampled  at  second  intervals.  It 
is  then  multiplied  with  the  complex  conjugate  of  the  CSI  represented  by  associ¬ 

ated  with  that  time  interval. 
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Combiner 


Figure  6.  BPSK  with  time  diversity  (After  Referenee  [9].)- 

Let  ,  where  A:  =  0,1 ,  be  the  transmitted  baseband  BPSK  signal  defined  as 

x„  =  x„  (t)  =  =  A  for  0  <  t  f 

.  (2.136) 

X,  =  X|  (t)  =  Ae  =-A  for  0  <  t  <  7). 

The  energy  £  of  a  transmitted  time  diversity  signal  is  [10] 

E  =  |xo  {t)f  dt  =  |xi it)f  dt  =  A% .  (2.137) 

As  the  transmitted  signal  passes  through  the  multipath  fading  ehannel  in  one  time  interval 
7)  seeonds,  it  is  multiplied  by  a  ehannel  eoeffieient  and  AWGN  (t)  is  added. 

Thus,  the  reeeived  diversity  signal  (t)  in  one  time  interval  is  given  by 

y,{t)  =  h,e^^'x,  for0<t<7;  ,  (2.138) 

where  /  is  the  time  diversity  reeeption  and  A:  =  0,1 . 

The  symbol  represents  the  signal  after  the  eorrelation  reeeiver.  The  signal  is 
given  by 

J  +  A,  when  A:  =  0  (2.139) 

\-A^T^hie^‘^'  +  A;  when  A:  =  1, 
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where  N,  is 


N,  =  n,  {t^dt . 


(2.140) 


Substituting  the  energy  E  defined  in  Equation  (2.137)  into  Equation  (2.139),  we  get 


-EhiC^'^'  +N, 


when  k  =  0 
when  k  =  l. 


(2.141) 


The  signal  after  multiplieation  with  CSl  of  the  /-th  reeeption  is  given  by 


Ehf+h,e-J^‘Ni 

-Ehf+hie-^'^'Ni 


when  A:  =  0 
when  A:  =  1 . 


(2.142) 


The  symbol  z^  represents  the  signal  eombined  with  all  L  reeeptions.  The  signal  z^.  is 


given  by 


+  h^e  when  ^  =  0 

i=\ 

+  whenA:  =  l. 


(2.143) 


The  eomparator  eompares  the  real  part  of  z^  to  zero.  If  Re[z^]  >  0  ,  a  binary  1  is 
the  output.  If  Re [z^.]  <  0,  a  binary  0  is  the  output.  Then  Re[z^.]  is  given  by 


Re[z,]  =  - 


Yj[Eh^  +Re[V^^'^'iV;]) 


when  A:  =  0 


+Re[^/i;e  when  A:  =  1. 


(2.144) 


The  signal  Re[z^]  is  a  Gaussian  random  variable  for  a  fixed  set  of  |/2,e  |  [10].  The 

probability  of  an  error  is  given  by  [9,10] 


Pt=Q 


(2.145) 


34 


where  is  the  positive  mean  of  the  random  variable  Re[z^] ,  i.e.,  the  expected  value  of 
Re[z^]  when  is  transmitted  and  is  the  variance  of  Re[z^] .  The  Q-function 
Q{x)  is  given  by  [5] 


(2.146) 


Thus,  to  find  /), ,  it  is  necessary  to  find  and  .  The  positive  mean  Z^  is  given  by 

r-'LEhj 

(2.147) 

=  Ethf. 

l=\ 

The  baseband  equivalent  noise  (t)  is  complex  sum  of  two  IID  Gaussian  noise 
sources  given  by 

ni{t)  =  x,{t)  +  jy,{t).  (2.148) 

The  autocorrelation  functions  of  the  baseband  equivalent  noise  and  its  real  and  imaginary 
parts  are  defined  as  [5] 

Kn{^)  =  E[n,{t)n,{t  +  T)] 

^x.(^)  =  e[x,  (t)x,(t  +  r)]  (2.149) 

The  power  spectral  density  of  the  baseband  equivalent  noise,  or  its  real  and  imaginary 
parts,  is  the  Fourier  transform  of  their  respective  autocorrelation  functions.  The  two- 
sided  baseband  equivalent  power  spectral  density  (/)  for  each  diversity  reception  is 
defined  as  [10] 

K  \f\  ^  - 

^.(/)=  I  (2-150) 

0  |f|>-. 
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where  B  is  an  arbitrary  bandwidth  that  is  B>\lT^ .  The  power  speetral  density  of  the 
corresponding  real  bandpass  noise  *S'(/)  centered  about  a  carrier  frequency  is  [10] 


^(/)  = 


4 

E, 

4 

0 


|/+/Jsf 

otherwise. 


(2.151) 


The  variance  of  the  baseband  equivalent  noise  cr^  is  equal  to  the  total  power  of 
the  noise  given  by  [5] 


o-:  =  NB . 


(2.152) 


The  total  power  of  the  baseband  equivalent  noise  is  the  sum  of  the  real  and  imaginary 
parts  [5].  Let  cr^  represent  the  variance  of  the  real  part  and  of  the  imaginary  part. 
Then  the  total  power  of  the  baseband  equivalent  noise  is  given  by  [5] 


^2  _  ^2 

nx/ 


(2.153) 


The  real  and  imaginary  parts  of  the  noise  are  IID  with  equal  variance.  Then  from  Equa¬ 
tion  (2.153)  the  variance  of  the  real  part  cr^  can  be  derived  by 


rx"  =  2ct" 

nx/ 

=NB 

nxi  o 

nx,  n. 


(2.154) 


Then  the  power  spectral  density  of  the  real  part  of  baseband  equivalent  noise  is  given  by 


^J/)  = 


2 


l/N 


B 


|/|>f 


(2.155) 


The  noise  power  after  the  matched  filter  and  multiplication  by  the  channel  parameter  is 
given  by  [7,10] 
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cr,  = 


Nht 


£|^o(/)f#. 


Using  Parsevel’s  theorem  in  Equation  (2.156),  we  get  [5] 

NX 


=- 


dt 


nXe 


(2.156) 


(2.157) 


The  varianee  is  the  sum  of  the  noise  power  in  eaeh  diversity  reeeption  and  is  given 
by  [9] 


cr, 


Re[z 


/=1 

L  \T  jX  ] 


=z 


NXE 


(2.158) 


/=1 


i=\ 


Using  Equations  (2.147)  and  (2.157)  in  Equation  (2.145),  we  obtain  [9,10] 

Ef^hf 


P,=Q 


=  Q 


l=\ 


\xx 

^  i=\ 

I— yy 


(2.159) 


Equation  (2.159)  is  a  funetion  of  a  random  variable  .  Eet  y  be  defined  as 


Then  Equation  (2.159)  beeomes 


E  1  2 

-'''o  1=1 


=  e(^/2r). 


(2.160) 


(2.161) 
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Note  that 


Yi  = 


E}^ 


(2.162) 


represents  the  instantaneous  SNR  in  the  /-th  reeeption  [10]. 

In  the  course  of  developing  the  multipath  model,  the  distribution  for  hj  is  already 
known.  If  the  random  variable  R)  is  defined  as 


W,=H^ 


(2.163) 


then  the  probability  distribution  function  is 


(2.164) 


Note  that  the  average  SNR  per  diversity  reception  is  given  by 


//=E[/,]=E 


Eh^ 


^  T7ri,2l  E  r  1  E2a^ 


iVn 


(2.165) 


The  characteristic  function  of  R)  is  given  by  [10] 


^w,  (®) 


1 

(l-72cr"®) 


(2.166) 


The  probability  distribution  function  of  the  sum  of  L  IID  random  variables  is  the 
T-fold  convolution  of  the  original  probability  distribution,  or  the  multiplication  of  L 
characteristic  functions  [5,10].  If  If  represents  the  sum  of  L  random  variables  W, ,  then 
the  characteristic  function  of  W  is 


Fw  (®)  = 


(2.167) 


This  is  the  characteristic  function  of  the  chi-square  distribution  of  degree  2L .  Therefore, 
the  probability  distribution  function  is  given  by 
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fw{w)=- 


w 


.i-1 


2V"^r(i) 

The  parameter  L  will  always  be  an  integer  so 

r(T)  =  (T-l)!. 

Then  Equation  (2.168)  becomes 


u{w) . 


Transforming  this  distribution  with 


2V"^(T-1)! 


f 

r  =—w, 


(2.168) 


(2.169) 


(2.170) 


(2.171) 


where 


dy  _  E 
dw  N„ 


(2.172) 


we  get 


/r(r)  = 


1  w" 


-w/ 2a 


^2V"^(T-1)! 

iVn 


(2.173) 


Equation  (2.173)  becomes 


fr{r)  = 


r 


L-\ 


y 

E2(J^ 


^E2a^ 


V  ^0  y 


(E-1)! 


(2.174) 


Substituting  the  average  SNR  y, ,  Equation  (2.165),  into  Equation  (2.174),  we  obtain 


fr{r)  =  — 


r 


.L-l 

-e^‘u 


ir)  (^-1)! 


(/)• 


(2.175) 
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Equation  (2.175)  is  the  probability  distribution  function  for  E  as  a  function  of  y . 
The  probability  of  bit  error  is  also  written  as  a  function  of  y .  The  average  probability  of 
bit  error  in  a  multipath  fading  channel  can  be  derived  by  taking  the  expectation  of  the  bit 
error  probability,  that  is  [9,10], 


P.  =  KPAr)] 


Substituting  Equations  (2.161)  and  (2.175)  into  (2.176),  we  get 


R  = 


r 


L-\ 


(r,)  (i-i)! 


ri 


dy . 


(2.176) 


(2.177) 


This  has  a  solution  given  by  [10] 

(l-n) 


P.  = 


where  the  parameter  u  is  given  by 


^(L-l  +  l^l  (1  +  m) 


1=0 


u  = 


/ 


(2.178) 


r 


'1  +  y 


(2.179) 


In  time-diversity,  the  combiner  adds  L  diversity  receptions  in  succession.  The 
transmit  interval  7)  is  related  to  the  bit  duration  7],  by  [9] 


n=LT^. 


(2.180) 


Eet  the  energy  per  bit  be  defined  as 

e,  =  aX. 


(2.181) 


The  correlation  filter  of  the  MRC  integrates  over  one  diversity  reception  7) .  The  energy 
per  bit  relates  to  energy  per  reception  by  [9] 

E,=LE,  (2.182) 


or 


E  = 


(2.183) 
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Substituting  Equation  (2.183)  into  Equation  (2.165),  we  get 


Yr 


LN, 


(2.184) 


In  spaee-diversity,  the  MRC  is  eonfigured  as  in  Eigure  7.  The  spaee-diversity  re- 
eeptions  eome  from  separate  reeeivers,  and  the  time  interval  T  of  the  matehed  filter  is  one 
bit  period  7), ,  i.e.,  T^=T  .  Therefore,  the  energy  per  diversity  reeeption  is  equal  to  the 
energy  per  bit 


E  =  E,. 

The  average  SNR  per  reeeption  beeomes 

-  E.lcr^ 

r,=— — 

'  ^  -\T 


(2.185) 


(2.186) 


Eigure  7.  MRC  spaee  diversity  eonfiguration  (After  Referenee  [9].). 
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Of  special  interest  is  when  L-\,  i.e.,  no  diversity.  The  average  SNR  per  recep¬ 
tion  for  both  time  and  space  diversity  schemes  reduces  to 


Yr 


(2.187) 


The  probability  of  bit  error  reduces  to 


P  = 


(1-^) 

2 


(2.188) 


Substituting  Equation  (2.179)  into  Equation  (2.188),  we  obtain 


f  —  ^ 

1 

1 

2 

1  >7 
-1- 

i 

(2.189) 


This  result  is  the  probability  of  bit  error  P^  for  a  BPSK  SISO  system  in  a  multipath  fad¬ 
ing  channel  [9,10].  Equation  (2.189)  will  be  the  benchmark  for  simulations  in  a  multi- 
path  fading  environment.  Eurthermore,  the  means  to  analyze  an  MRC  receiver  is  ob¬ 
tained.  This  will  be  essential  in  the  analysis  of  the  design. 


F.  SPACE-TIME  CODING 

“An  effective  and  practical  way  to  approaching  the  capacity  of  MIMO  wireless 
channels  is  to  employ  space-time  coding”  [1].  Space-time  coding  is  simply  transmitting 
signals  over  the  channel  from  multiple  antennas  to  multiple  antennas  in  such  a  way  that 
capacity  is  maximized.  Maximum  capacity  can  only  be  approached.  Space-time  codes 
are  designed  to  achieve  both  time  and  space  diversity  to  improve  the  error  performance 
and  achieve  diversity  gain  [1]. 

Another  mechanism  to  improve  performance  is  error  control  coding  (ECC)  [1]. 
Error  control  coding  arranges  the  transmitted  symbols  in  such  a  way  as  to  increase  im¬ 
munity  to  noise.  It  is  realized  by  introducing  redundant  bits  to  the  information  bit  stream, 
allowing  the  receiver  to  detect  and  possibly  correct  errors.  However,  increasing  the 
number  of  bits  in  the  bit  stream  increases  the  bandwidth  required  to  transmit  the  informa¬ 
tion  [10,1 1].  If  the  error  control  coding  is  designed  in  conjunction  with  the  modulation 
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technique  and  the  transmit  diversity  scheme,  then  it  is  possible  that  there  is  no  bandwidth 
expansion.  With  ECC,  coding  gain  is  achieved  [1]. 

Overall,  it  is  possible  to  achieve  both  diversity  and  coding  gain  with  space-time 
codes  [1].  (The  general  effects  of  diversity  and  coding  on  performance  are  illustrated  in 
Figure  8  from  Reference  [13].).  In  general  it  depicts  the  possibility  of  improving  per¬ 
formance  by  the  introduction  of  diversity  and  ECC.  The  focus  of  this  thesis  is  diversity 
gain;  the  addition  of  error  control  coding  is  left  for  future  work. 


Figure  8.  Diversity  and  coding  gain  (From  Reference  [13].). 

In  this  chapter,  the  fundamentals  necessary  to  understand  a  MIMO  system  in  a 
multipath  fading  environment  were  introduced.  It  was  seen  that  the  capacity  of  a  MIMO 
system  can  increase  linearly  with  the  number  of  transmit  and  receive  antennas.  To  ap¬ 
proach  this  capacity,  space-time  codes  are  designed  to  incorporate  diversity  and  error 
control  techniques.  A  baseband  model  of  a  slow  fading,  frequency  nonselective  multi- 
path  fading  channel  was  obtained  for  simulation.  Easily,  diversity  techniques  with  the 
MRC  receiver  were  examined  to  facilitate  analysis  of  the  design.  In  the  next  chapter,  the 
Alamouti  space-time  coding  scheme  is  introduced,  and  a  transmitter  and  receiver  design 
will  be  developed. 
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III.  MULTIPLE-INPUT-MULTIPLE-OUTPUT  TRANSMITTER 
AND  RECEIVER  USING  THE  ALAMOUTI  SPACE-TIME  CODING 

SCHEME 


This  chapter  introduces  the  Alamouti  space-time  code.  Once  the  code  is  pre¬ 
sented,  it  will  be  used  to  design  the  MIMO  transmitter  and  receiver  using  BPSK  modula¬ 
tion.  The  designs  will  be  simulated  in  Chapter  IV  and  bit  error  performance  will  be  ana¬ 
lyzed  and  tested. 


A,  THE  ALAMOUTI  SPACE-TIME  CODING  SCHEME 

The  Alamouti  space-time  scheme  was  introduced  in  Reference  [14].  It  is  a 
space-time  code  for  transmission  through  two  antennas,  N^=2.  It  will  be  used  in  a 

MIMO  configuration  with  two  receive  antennas,  =  2. 

To  describe  signals  more  precisely  from  transmitter  to  receiver,  the  following  no¬ 
tation  will  be  adhered  to  through  the  remainder  of  this  thesis.  Let: 

•  x"  represents  the  signal  x  in  the  n-th  transmission  interval  from  the  m-th 
antenna, 

•  x"'  represents  the  signal  over  the  complete  code  sequence  from  the  m  -th 
antenna, 

•  square  brackets  [  ]  encompass  the  total  space-time  code.  The  elements 

within  the  brackets  separated  by  commas  represent  the  code  signal  of  one 
time  interval,  and 

•  asterisk  *  represents  the  complex  conjugate  of  the  signal. 

Let  X]  and  Xj  represent  the  baseband  equivalent  of  signals  for  BPSK  modulated 
bits  over  two  consecutive  transmission  intervals  of  T  seconds.  Each  is  represented  by 
±A  .  Let  x'  and  be  the  Alamouti  space-time  coded  signals  transmitted  from  antenna 
one  and  antenna  two,  respectively.  Then  the  Alamouti  space-time  coded  signals  are  de¬ 
fined  by  [1,13] 
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(3.1) 


X 


X 


=[xi,-x;] 

=[x2,x;]. 


Although  Xj  and  Xj  are  real  in  this  ease,  the  eomplex  eonjugate  remains  for  generality 
and  possible  transformation  to  eomplex  values  in  future  work. 

The  Alamouti  spaee-time  eode,  as  it  emanates  from  either  antenna,  eontains  in¬ 
formation  neeessary  to  demodulate  sueeessive  information  bits.  The  eombination  of  the 
two  eonseeutive  signals  from  both  transmit  antennas  eorresponds  to  time  diversity.  Two 
antennas  are  used  at  the  reeeiver,  thus,  spaee  diversity  is  also  aehieved  [1]. 


The  Alamouti  scheme  has  two  important  features  [1].  First,  the  coded  signals 
from  the  two  transmit  antennas  are  orthogonal  over  the  code  interval.  If  x'  and  x^  are 
considered  as  complex  vectors  then  they  are  orthogonal  if  their  inner  product  equals  zero 
[4,6].  Specifically, 


Expansion  of  left  hand  side  of  Equation  (3.2)  gives 


1  (  2  r  * 
X  I  X  1  =  Xj  -Xj 


=  XjXj  -  XjXj  =  0. 


(3.2) 


(3.3) 


Thus,  x'  and  x^  are  orthogonal. 


To  describe  the  second  feature,  first  let  the  code  matrix  be  defined  by  X .  The 
code  matrix  X  is  the  combination  of  x'  and  x^  and  is  given  by  [1] 


Xi  -x’ 

— 

* 

X 

Xj  Xj 

(3.4) 
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Then,  XX*^  is 


XX’" 


XjXj  +  XjXj  XjXj  -  XjXj 

_X2X;-XiX2  X^x’+Xi’Xi 

“l  |2  I  |2 

|Xi  I  +  jXj  I  0 

0  1-^11^+ 1-^21^ 

(Kr+Kr)i2- 


(3.5) 


Equation  (3.5)  represents  a  scalar  multiplied  by  an  identity  matrix.  Note  as  well  that  the 
orthogonal  feature  of  the  Alamouti  code  is  encompassed  in  this  feature.  Furthermore, 
this  feature  represents  a  condition  similar  to  the  example  of  the  capacity  of  MIMO  sys¬ 
tems  in  Chapter  II.  It  was  found  that  with  this  matrix,  i.e.,  a  scalar  multiplied  by  the  iden¬ 
tity  matrix,  the  capacity  increased  linearly  with  the  rank  of  H .  This  feature  and  its  im¬ 
portance  will  become  clear  in  the  next  chapter. 

With  the  Alamouti  space-time  coding  scheme  defined,  the  transmitter  and  re¬ 
ceiver  can  now  be  designed.  A  generic  transmitter  and  receiver  is  described  in  Reference 
[1].  The  transmitter  and  receiver  described  next  is  the  implementation  using  BPSK 
modulation. 


B,  TRANSMITER 

The  objective  of  the  transmitter  is  clear;  produce  two  BPSK  signals  from  an  in¬ 
formation  bit  stream  for  transmission  from  two  transmit  antennas  according  to  Equation 
(3.1).  The  transmitter  design  is  shown  in  Figure  9.  Eet  the  binary  information  source  be 
represented  by  b  .  The  binary  information  source  undergoes  serial-to-parallel  conversion 
to  produce  two  binary  bits  and  b^  in  successive  intervals.  Next,  the  two  binary  sym¬ 
bols  are  modulated  by  baseband  BPSK  modulators  to  give  the  corresponding  two  base¬ 
band  equivalent  signals  Xj  and  Xj .  In  turn,  Xj ,  Xj  and  their  complex  conjugates  are 

routed  to  two  parallel-to-serial  converters  to  produce  x'  and  as  defined  by  Equation 
(3.1). 
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An  important  feature  of  the  transmitter  is  that  the  transmission  rate  of  the  infor¬ 
mation  source  is  equal  to  the  transmission  rate  of  each  of  the  two  transmitted  signals 

jc'  and  .  Each  bit  period  is  T  seconds  and  is  equal  to  the  original  bit  period.  The  se- 
rial-to-parallel  and  parallel-to-serial  conversion  combination  maintains  the  transmission 
rate  of  the  binary  information  source  .  This  is  another  important  feature  of  space-time 

coding.  In  MIMO  systems,  time  and  space  diversity  can  be  achieved  without  bandwidth 
expansion[l].  In  this  case,  the  Alamouti  space-time  code  achieves  both  space  and  time 
diversity  while  maintaining  the  original  transmission  rate  through  the  channel. 


Figure  9.  MIMO  Transmitter. 


C.  RECEIVER 

The  transmitter  was  relatively  straightforward  to  design.  The  receiver  is  some¬ 
what  more  complicated.  Let  the  signal  available  at  each  receive  antenna  be  the  repre¬ 
sented  by  r'  and  .  In  the  absence  of  AWGN  and  multipath  fading  (i.e.,  all  channel 
coefficients  equal  one),  r'  and  are  given  by 


r 

r 


1 

2 


=  jc'  + 


=  jc'  +JC^ 


(3.6) 


Incorporation  of  the  definitions  for  x'  and  x^  from  Equation  (3.1)  into  Equation  (3.6) 
gives 
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(3.7) 


r'=[x„-x;]  +  [x„x;] 
r^=[x„-x;]  +  [x2,x:]. 

Over  the  eode  sequence  the  same  information  is  available  at  each  receive  antenna.  This 
is  possible  because  of  the  assumption  that  all  channel  coefficients  are  equal  to  one.  In  a 
fading  environment  this  will  not  be  the  case.  Nevertheless,  there  is  time  and  space  diver¬ 
sity  capability,  an  advantage  that  should  be  exploited.  The  logic  of  the  receiver  design  is 
to  exploit  this  and  maximize  the  diversity  reception  of  Xj  and  Xj . 


Let  the  signals  be  isolated  in  one  transmission  interval  T.  A  serial-to-parallel 
converter  performs  this  function.  The  received  signals  by  transmission  interval  are 


r'  =[r;,riy[x,+x,,xl-xl]. 

(3.8) 

and 

1  V  2  2~\  V  *  *n 

(3.9) 

If  the  complex 

conjugate  is  taken  on  and  ,  the  received  signals  become 

r/  =Xi+X2 

[rl]  =Xi-X2, 

(3.10) 

and 

q"  =Xi+X2 
[r^)  =Xi-X2. 

(3.11) 

From  a  closer  look  at  Equations  (3.10)  and  (3.11),  there  is  a  total  of  four  signals 
available.  The  signals  in  each  time  interval,  distinguished  with  different  subscripts,  rep¬ 
resent  time  diversity  signals.  The  signals  received  by  each  antenna,  distinguished  by  dif¬ 
ferent  superscripts,  represent  space  diversity  signals.  By  combining  correctly,  both  Xj 

and  Xj  can  be  isolated.  Let  and  z"’ ,  where  m  =  1, 2  ,  be  the  time  diversity  combina¬ 
tion  on  the  m-th  space  diversity  reception  that  isolates  Xj  and  Xj  respectively.  Then  the 
time  diversity  combination  for  Xj  is  given  by 
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(3.12) 


Lastly,  the  signals  Zj  and  Zj  need  only  be  demodulated  to  produee  and  , 
where  the  distinguishes  the  demodulated  binary  signal  from  the  original  signal. 

Lastly,  hj  and  4  are  sent  to  a  parallel-to-serial  converter  to  reproduce  the  demodulated 
information  bit  stream  b  .  The  receiver  design  is  illustrated  in  Figure  10. 

The  MIMO  system  that  features  the  Alamouti  space-time  code  is  now  designed. 
The  system  incorporates  transmit  and  space  diversity  techniques  in  order  to  combat  the 
effects  of  the  multipath  fading  channel.  In  the  next  chapter,  the  design  is  implemented  in 
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Matlab  with  Simulink  and  is  tested  for  performance  with  AWGN  alone  and  AWGN  in  a 
multipath  fading  environment. 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


52 


IV.  ANALYSIS  AND  SIMULATION 


In  this  chapter  the  error  performanee  analysis  and  simulation  of  the  MIMO  system 
with  the  Alamouti  spaee-time  code  is  presented.  The  system,  as  deseribed  in  Chapter  III, 
was  simulated  in  Matlab  using  Simulink.  The  system  is  analyzed  and  simulated  in  a  pro¬ 
gressive  manner  from  an  AWGN  channel  to  the  multipath  fading  ehannel.  To  faeilitate  a 
smooth  transition  from  each  progressive  step,  the  changes  to  the  original  simulation  setup 
are  highlighted.  The  simulation  of  the  transmitter  and  reeeiver  are  described  first. 

For  ease  of  presentation,  the  Matlab  simulation  bloeks,  as  they  are  used  in  the 
simulations,  are  presented  in  Appendix  B.  Further  information  eoncerning  the  bloeks  and 
the  parameters  are  found  in  Referenee  [15].  The  parameters,  as  depleted  in  Appendix  B, 
are  maintained  throughout  all  simulations. 

Each  simulation  is  compared  to  a  baseband  equivalent  SISO  BPSK  system.  An 
important  parameter  for  eomparison  is  the  average  energy  per  bit  .  In  a  single  ehannel 

BPSK  system,  is  equal  if  a  binary  bit  zero  or  one  is  transmitted.  Let  the  modulated 
signal  for  a  binary  zero  be  represented  by  Xg  (t) .  If  the  amplitude  over  a  period  of  T 
seeonds  is  represented  by  A  ,  then  is  given  by  [10] 

E,=\l\x,{t)\dt  =  A^T.  (4.1) 

The  energy  per  bit  is  referenced  to  the  output  of  the  transmitter,  thus  the  transmitted 
power  into  the  ehannel  for  the  SISO  system  is  defined  by 

(4.2) 

The  symbol  is  defined  as  the  energy  per  symbol  and  is  used  to  represent  the  energy 
per  symbol  period  of  the  MIMO  system.  This  is  done  in  order  to  distinguish  between  the 
energies  of  the  SISO  and  MIMO  systems.  The  mathematical  representation  of  E^  will  be 

derived  later  in  this  ehapter.  It  will  also  be  refereneed  to  the  output  of  the  transmitter, 
and  the  transmitted  power  into  the  channel  for  one  MIMO  antenna  is  defined  by 
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(4.3) 


p 

^  MIMO  j,  ■ 

A,  SIMULATION  OF  THE  MIMO  TRANSMITTER  AND  RECEIVER 

1.  MIMO  Transmitter  Simulation 

The  block  diagram  of  the  transmitter  simulation  is  given  in  Figure  1 1 .  The  trans¬ 
mitter  information  source  is  a  Bernoulli  binary  generator  block.  It  produces  the  binary 
information  source  b  .  The  output  is  a  binary  bit  stream  of  Is  and  Os  according  to  three 
parameters.  The  first  is  the  probability  of  zero.  Since  the  ones  and  zeros  are  equally 
likely,  this  parameter  is  set  to  0.5.  The  second  parameter  is  the  initial  seed  for  the  ran¬ 
dom  number  generator.  This  value  is  a  prime  number  that  must  be  different  from  all 
other  random  sources  in  the  simulation.  By  ensuring  different  seeds,  the  statistical  inde¬ 
pendence  of  all  sources  is  simulated  [15].  The  last  parameter  is  sample  time,  which  is  set 
at  the  value  T  [15]. 


r-LJv_iT_rv 

Bernoulli 

b 

ll 

-♦ 

Ciieinterlacpr 

Binary 

Bernoulli  Binai 
Generator 


OOCLZ 

1 

DrOfN 

* 

BPSK  1 

j 

[A. 

Modulator 

^1 

1  1 

Baseband  1 

Math 

“v_rWvw 

BPSK 


BPSK 

Modulator 

Baseband2 


FuncDonl 


■HZhG> 

Math  Gam 
Funcbon2 


Interlacer 


x* 


U 

Out) 


Interlacer  1  Unbutferl 


Interlacer 


Ineerlacer2  UnOufferZ 


Out2 


Figure  11.  MIMO  Transmitter  Simulation. 

The  next  two  blocks  in  the  MIMO  transmitter  are  the  buffer  and  deinterlacer.  The 
buffer  takes  the  input  bit  sequence  b  and  produces  a  2-by-l  frame  at  a  sample  time  of 
2T  seconds.  The  deinterlacer  splits  the  frame  into  two  separate  outputs,  isolating  b^  and 


54 


&2  for  manipulation.  Together,  the  buffer  and  deinterlaeer  provide  the  serial-to-parallel 
conversion  of  the  bit  stream.  The  buffer  must  wait  for  2T  seconds  before  producing  an 
output  that  was  generated  by  the  Bernoulli  binary  generator.  This  causes  a  delay  of 
2T  seconds. 

Now  that  hj  and  hj  are  separated,  each  is  sent  to  a  baseband  equivalent  BPSK 
modulator  to  produce  the  modulated  symbols  Xj  and  Xj .  The  modulation  of  the  binary 
data  is  accomplished  according  to  Table  1. 


Input  hj  2 

Output  Xj2 

0 

A 

1 

-A 

Table  1.  Modulation  of  binary  data. 


where  the  amplitude  A  =  1 .  The  modulated  signal  Xj  is  routed  through  a  complex  conju¬ 
gate  block  to  produce  x* .  Similarly,  Xj  is  routed  through  a  complex  conjugate  block  and 
a  negative  gain  block  to  produce  -x^ .  The  baseband  modulated  symbols  Xj ,  x* ,  Xj  and 
-Xj  are  then  routed  to  the  input  of  two  interlacer  blocks  to  form  the  final  outputs.  The 

interlacer  and  unbuffer  blocks  provide  the  parallel-to-serial  conversion.  The  output  of 
each  unbuffer  block  has  a  period  of  T  seconds.  The  output  identified  in  Figure  1 1  as 
“Outl”  is  the  code  sequence  for  antenna  1,  x' .  Likewise,  the  “Out2”  is  the  code  se¬ 
quence  for  antenna  2,  .  Each  conforms  to  the  required  code  sequence  of  the  Alamouti 

space-time  code  given  in  Chapter  III. 

2.  MIMO  Receiver  Simulation 

The  block  diagram  of  the  receiver  simulation  is  shown  in  Figure  12.  The  inputs 
indicated  in  the  figure  as  “Ini”  and  ‘Tn2”  represent  the  received  signals  from  antennas  1 
and  2,  respectively,  that  is  r'  and  .  The  buffer  and  deinterlacer  blocks  are  identical  to 
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the  transmitter  and  serve  as  the  serial-to-parallel  eonversion  of  the  reeeived  signals  r' 
and  .  However,  it  is  important  to  note,  as  with  the  transmitter,  the  buffer  bloeks  intro¬ 
duce  another  delay  of  2T .  The  total  delay  through  the  system  is  now  4T .  The  math 
function  blocks  implement  the  complex  conjugate  of  their  respective  inputs.  The  sum¬ 
ming  blocks  are  the  time  and  space  diversity  combining  described  in  Chapter  III. 


Function2 


Figure  12.  MIMO  Receiver  Simulation. 


The  baseband  equivalent  demodulator  blocks  demodulate  the  incoming  signals 
according  to  Table  2  [15]. 


Input  Zj2 

Output  hj  2 

Re[zi_2]>0 

0 

Re[zi_2]<0 

I 

Table  2.  Demodulation  to  binary  data 
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The  output  of  the  two  demodulator  blocks  are  the  binary  bit  streams  and  .  The  final 
two  blocks,  the  interlacer  and  unbuffer,  serve  as  the  parallel-to-serial  converter  producing 
the  demodulated  binary  information  b  . 

B,  SIMULATION  IN  AWGN 

The  MIMO  system  is  first  subjected  to  AWGN  only.  In  this  simulation,  it  is  as¬ 
sumed  there  is  no  fading.  This  implies  that  the  channel  coefficients  h^j  are  all  constant 

and  equal  to  one.  A  block  diagram  of  the  MIMO  system  in  AWGN  is  illustrated  in 
Figure  13.  Both  the  MIMO  transmitter  and  receiver  are  collapsed  into  subsystems  for 
ease  of  presentation.  They  remain  the  same  as  those  described  in  Figure  1 1  and  Figure 
12,  respectively.  In  this  simulation,  two  gain  blocks  are  introduced  at  the  outputs  of  the 
transmitter.  This  is  done  to  normalize  the  total  power  transmitted  into  the  channel  for 
comparison  to  a  SISO  BPSK  system.  The  MIMO  channel  in  AWGN  consists  of  the  out¬ 
put  of  the  two  gain  blocks  to  the  two  inputs  of  the  MIMO  receiver.  The  two  AWGN 
channel  blocks  are  complex  white  Gaussian  noise  sources  that  add  to  the  incoming  sig¬ 
nals.  They  are  represented  by  and  . 


Channel2 


F igure  13.  S  imulation  in  AW GN . 
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1,  Simulation  Parameters 

The  factor  K  in  the  gain  blocks  is  determined  from  the  discussion  of  the  MIMO 
model  in  Chapter  II.  The  power  emanating  from  both  of  the  transmit  antennas  is  equal. 
Let  the  power  from  any  antenna  be  represented  by  .  Let  the  total  power  transmitted 
into  the  channel  be  represented  by  P  .  The  power  from  any  antenna  is  the  total  power  di¬ 
vided  by  the  number  of  transmit  antennas,  in  this  case  =  2.  The  power  from  any  an¬ 
tenna  P^"‘  is  given  by 


P 


Ant 


(4.4) 


The  energy  in  any  symbol  Xj ,  Xj  x*  or  xj  is  equal  whether  a  binary  zero  or 

one  is  transmitted.  The  distinction  between  the  SISO  and  MIMO  system  is  the  MIMO 
system  symbols  include  the  multiplication  of  K  by  the  gain  blocks.  Once  again,  let  the 
modulated  signal  for  a  binary  zero  be  represented  by  Xq  (t) .  If  the  amplitude  over  a  pe¬ 
riod  of  T  seconds  is  represented  by  A  ,  then  is  given  by  [10] 


E^  =  |^^|Xxo  {t)^dt  =  f^K^A^dt  =  K^A^T . 

The  power  emitted  by  each  antenna  is  equal  to  the  symbol  power  given  by  [7] 

p 

pAnt  _  ^ s  _  ^2^2 

~  T  ~ 


(4.5) 


(4.6) 


When  Equations  (4.4)  and  (4.6)  are  equated,  the  total  signal  power  into  the  channel  as  a 
function  of  amplitude  and  the  gain  factor  K  is  given  as 


—  =  K^A^ 

2 

P  =  2K^A\ 


(4.7) 


The  total  power  at  the  input  of  either  of  the  AWGN  blocks  with  no  fading  is  also 
given  by  P  .  This  is  set  to  1.0  Watt.  The  amplitude  A  is  1.0  Volt.  Using  this  informa¬ 
tion  in  Equation  (4.7),  we  obtain 


58 


(4.8) 


-- 

K  = 


Thus,  the  gain  factor  K  is  set  to  1/V2  in  this  and  all  subsequent  simulations.  The  energy 
per  symbol  of  the  MIMO  system  is 


E^  = 


A^T 

~Y~ 


(4.9) 


The  energy  per  symbol  for  Xj ,  Xj ,  x’  or  Xj  is  given  by  Equation  (4.9).  It  incorporates 
the  gain  factor  K .  Thus,  if  a  bit  zero  is  transmitted  in  the  MIMO  system,  represented  by 
Xq  (t) ,  it  is  given  by 


A 


Xg  (t)  =  for  0  <  t  <  r  . 

y/2 


(4.10) 


Comparing  the  energy  per  bit  of  the  SISO  system.  Equation  (4.1),  and  the  energy  per 
symbol  E^  of  the  MIMO  system.  Equation  (4.9),  we  have 


The  AWGN  channel  blocks  require  several  parameters:  an  initial  seed,  the  energy 
per  symbol-to-noise  ratio  EjN^  ,  the  power  of  the  input  signal,  and  the  symbol  period. 
The  symbol  period  is  set  to  T  seconds.  The  average  energy  per  symbol-to-noise  ratio 
EjNg  increases  with  each  simulation  run.  The  initial  seed  is  a  prime  number  that  is  dif¬ 
ferent  in  each  block  as  well  as  different  from  the  initial  seed  of  the  Bernoulli  binary  gen¬ 
erator  block. 

Einally,  the  error  rate  calculation  block  compares  each  bit  from  the  original  bit 
stream  b  with  the  demodulated  bit  stream  b  over  the  simulation  run.  It  calculates  the  bit 
error  rate  or  simulated  probability  of  bit  error.  The  total  delay  of  AT  is  incorporated  in 
this  block  to  ensure  that  the  comparison  of  b  and  b  is  synchronized  in  time. 
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2,  Performance  Analysis 

At  each  receive  antenna,  the  received  signals  r'  and  are  given  by 

r'  =(x'  +  + 

=(x'  +x^  j  + 


(4.12) 


The  complex  valued  vectors  «'  and  in  Equation  (4.12)  represent  the  noise  over  two 
consecutive  bits  added  at  receive  antennas  1  and  2,  respectively.  If  the  definitions  for  x' 
and  x^  from  Chapter  III  are  incorporated,  then  Equation  (4.12)  becomes 


r'  =  ([xi , -x* ]  +  [xj, Xi  ])  +  [n\ , n\ ] 

=  ([-^1 ,  --^2  ]  +  [-^2  Ai*  ])  +  \nl ,  «2  ]  • 


(4.13) 


After  serial-to-parallel  conversion,  the  signals  in  each  time  and  space  diversity  reception 
are  given  by 


and 


iriinf/  \  \  (  *  *\ 

r  =[ri  ,r2j  =  [(Xi+X2)  +  ni,(Xi -Xjj 


~  (-^1  +X2)  +  nl,{xl  -x^'^  +  nl 


(4.14) 


(4.15) 


The  time  diversity  combination  of  the  received  signal  on  antenna  1 ,  represented 
by  z\  and  Zj ,  respectively,  are 


and 


+(fi) 

=  {x,+x^)  +  n\  +((xi -x’j  +  n’) 
=  2xj  +n\  +^2*’ 


=  (xi +X2)  +  n! -((x* -x’j  +  n’) 

=  2x2  +^1  ■ 


(4.16) 


(4.17) 
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Similarly,  the  time  diversity  eombination  of  the  received  signal  on  antenna  2,  zf  and  z\  , 
are 


2  ,2,2* 
Zj  =  2Xj  +  +  ^2 

Zj  =  2^2  +  «!  -  ^2  • 


(4.18) 


The  space  diversity  combining  of  Xj  and  X2  are  represented  by  Zj  and  Z2 .  After 
space  diversity  combining,  Zj  is 


Zj  =  z\  +  zf 

=  2xj  +  n\  +  +  nf  +  (4.19) 

A  1  1*  2  2* 

and  Z2  is 

1  ,  2 

^2  =  Z2  +  Z2 

=  2X2  +  n\  -  +  2X2  +  nf  -  (4.20) 

A  1  1*  2  2* 

=  4X2  +  -  ^2  +  ^1  “  ^2  • 

It  is  important  to  note  the  symmetry  of  the  receiver.  In  Equations  (4.19)  and 
(4.20),  Xj  and  X2  are  the  BPSK  baseband  modulation  of  two  binary  random  variables  that 

are  IID.  The  noise  components  and  represent  zero-mean  IID  complex 

Gaussian  random  variables.  The  sum  of  the  means  is  zero,  and  the  sum  of  variances  is 
equal  in  both  Equations  (4.19)  and  (4.20)  [5].  Thus,  the  probability  of  error  in  demodu¬ 
lating  Zj  and  Z2  are  equal.  The  overall  probability  of  bit  error  of  the  system  is  the  aver¬ 
age  of  the  two.  Since  they  are  equal,  the  probability  of  bit  error  of  the  system  is  the  prob¬ 
ability  of  error  in  demodulating  either  Zj  or  Z2 . 


Eet  the  overall  probability  of  bit  error  be  represented  by  .  To  determine  ,  Zj 

is  used.  The  demodulator  takes  the  real  part  of  Zj .  Eet  z  be  the  real  part  of  Zj ,  given  by 

z  =  Re(zj )  =  Re(^4xj  +  n\  +  nf  +  nf  +  ri^  j 

=  4xj  -l-Re(nJ^-l-Re(^n2*)  +  R^(^f  )“R^(^r)- 
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Let  the  sum  of  the  real  parts  of  all  the  noise  eomponents  be  represented  by  n  .  Then 
Equation  (4.21)  reduees  to 

z  =  4xj+n.  (4.22) 

Equation  (4.22)  defines  a  Gaussian  random  variable  Z  .  The  probability  of  bit  error  for  a 
Gaussian  random  variable  [9,10]  is 


P.  =  Q 


f 


(4.23) 


where  is  the  positive  mean  of  Z  (i.e.,  the  mean  of  Z  when  a  binary  0  is  transmitted). 
Eet  Xq  (t)  represent  the  modulated  signal  Xj  when  a  binary  0  is  transmitted.  If  matched 

filter  conditions  are  assumed,  then  the  positive  mean  Z^  at  the  output  of  the  receiver  is 
given  by  [9,10] 

=  A\lxl{t)dt  (4.24) 

=  4£.- 

The  power  spectral  density  of  the  real  part  of  each  noise  component  is 
Watts  per  Hertz  [15].  Eet  the  total  noise  power  of  the  real  part  of  the  noise  be  repre¬ 
sented  by  <jI  .  The  noise  power  cr^  for  each  noise  component  after  the  matched  filter  is 


cr„  = 


By  Parsevel’s  theorem  [5],  Equation  (4.25)  converts  to 


Jo  2 


(4.25) 


(4.26) 


Eet  the  total  noise  power  of  Z  be  represented  by  cr^ .  Then  cr^  is  sum  of  each  of  the 
noise  components  crl  and  is  given  by 
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(4.27) 


cr,  =  4cr„  =  NE^ . 


O  OS 


^  4E  ^ 


Using  Equations  (4.24)  and  (4.27)  in  Equation  (4.23),  we  get 

Pe=Q 

f 

=  Q 


4Ee, 


16E, 


^  J 
\ 


(4.28) 


vv  j 


Substituting  the  energy  per  bit  defined  in  Equation  (4.1 1)  into  Equation  (4.28),  we  obtain 

/ 

Pe=Q 


(4.29) 


=  Q 


2/2^. 


N, 


The  probability  of  bit  error  /),  for  a  baseband  equivalent  SISO  BPSK  system  is 
given  by  [9,10] 


P,  =  Q 


v^zy 


(4.30) 


In  this  ease  is  given  by 


and 


=  E. 


h  ’ 


(4.31) 


(4.32) 


Thus,  E  is 


P,  =  Q 


2E. 


(4.33) 


A  comparison  of  Equations  (4.29)  and  (4.33)  reveals  an  improvement  in  bit  error 
performance  in  the  MIMO  system  by  a  factor  of  two  over  the  SISO  system,  for  equal 
performance,  this  equates  to  a  3-dB  gain. 
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3,  Performance  Simulation  Results 

The  simulation  was  conducted  with  increasing  EjN^  .  The  period  T  was  set  to 

0.0001  seconds.  The  simulation  run  time  was  set  at  500.0  seconds.  Thus,  the  total  num¬ 
ber  of  bits  in  a  single  simulation  run  was  five  million  bits. 

The  simulated  probability  of  bit  error  for  the  MIMO  system  is  plotted  in  Figure 
14.  The  results  are  plotted  with  the  theoretical  results  obtained  in  Equation  (4.29).  In 
addition,  the  theoretical  probability  of  bit  error  for  a  baseband  equivalent  SISO  system  is 
plotted  for  comparison(Equation  (4.33)).  The  simulated  results  follow  the  theoretical  re¬ 
sults  very  well. 


Figure  14.  Results  in  AWGN. 

They  deviate  starting  at  approximately  5.0  dB  and  above,  but  only  by  approxi¬ 
mately  0.25  dB.  This  is  a  result  of  the  number  of  bits  in  a  single  simulation  run.  If  more 
bits  are  sent  for  comparison  in  a  single  simulation  run,  then  comparison  of  simulated  re¬ 
sults  with  the  theoretical  results  improves.  This  was  confirmed  with  several  experiments. 
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Nevertheless,  overall  the  simulation  of  the  MIMO  system  with  five  million  bits  follows 
the  theoretieal  analysis  of  bit  error  performanee  very  well. 


C.  SIMULATION  IN  A  MULITPATH  FADING  CHANNEL 

Next,  the  effeet  of  the  multipath  fading  ehannels  are  added.  The  simulation  bloek 
diagram  is  given  in  Figure  15.  A  multipath  fading  ehannel  is  inserted  for  eaeh  independ¬ 
ent  path  from  transmitter  to  reeeive  antenna.  There  are  a  total  of  four  independent  paths, 
eaeh  with  a  ehannel  eoeffieient  represented  by  h^J . 


1,  Simulation  Parameters 

The  multipath  fading  model  was  derived  in  Chapter  II.  Eaeh  multipath  fading 
model  is  identieal  with  respeet  to  eonfiguration.  A  bloek  diagram  is  given  in  Figure  16. 
The  model  is  the  eomplex  sum  of  two  IID  Gaussian  random  noise  sourees  with  zero- 
mean  and  equal  varianee.  To  distinguish  these  Gaussian  random  noise  sources  from  the 
AWGN,  they  will  be  called  Gaussian  channel  fade  noise  sources.  The  three  parameters 
of  the  Gaussian  channel  fade  noise  sources  that  must  be  determined  and  controlled  are  the 
variance,  the  sample  time  and  the  initial  seed.  The  initial  seed  determines  the  independ¬ 
ence  of  the  channel.  Each  of  the  seeds  within  the  model  is  different.  In  this  thesis  the  fo¬ 
cus  is  on  the  case  of  all  channel  coefficients  being  statistically  independent.  Therefore, 
each  of  the  initial  seeds  of  all  the  other  models  is  different.  Consequently,  there  are  eight 
different  seeds  for  the  Gaussian  channel  fade  noise  sources. 

Next  is  the  sample  time.  The  sample  time  is  determined  such  that  the  channel  co¬ 
efficients  remain  constant  over  the  code  sequence.  Thus,  the  channel  is  slowly  fading 
over  the  code  sequence.  The  duration  of  the  code  sequence  is  2T  seconds,  thus  the  sam¬ 
ple  time  is  set  to  2T  seconds. 
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Figure  15.  MIMO  in  multipath. 
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Figure  16.  Fading  model. 
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Lastly,  the  variance  of  the  Gaussian  channel  fade  noise  sources  is  set  to  a  value 
such  that  the  total  average  signal  power  entering  both  AWGN  blocks  is  1 .0  Watt.  To  de¬ 
termine  the  variance  of  the  noise  source,  the  energy  per  symbol  at  the  output  of  Multi- 
path  1  is  examined.  Let  the  energy  per  symbol  at  the  output  of  the  fading  channel  be  rep¬ 
resented  by  given  by  [9,10] 


M 

f  A^T 
~1 


dt 


(4.34) 


The  channel  coefficient  is  a  random  variable,  so  the  average  energy  per  symbol 
after  the  multipath  channel  is  the  mean  of  .  This  can  be  obtained  by  taking  the  ex¬ 
pected  value  of  Equation  (4.34)  [5].  Then  E^  is 


/2ii  A^T 


A^T 


(4.35) 


The  expected  value  of  is  the  second  moment  of  the  Rayleigh  distribution. 

The  second  moment  can  be  determined  from  the  mean  and  variance  of  the  Rayleigh  dis¬ 
tribution  defined  in  Chapter  II  by 


E[hl,~\  =  CTl+h^ 


( 

l—\ 

\n 

2 - 

+ 

1  2; 

[  ^ 

hj 

=  1ct\ 


(4.36) 


Note  that  <j^  represents  the  variance  of  the  Gaussian  channel  fade  noise  sources.  Using 
Equation  (4.36)  in  Equation  (4.35),  we  obtain  the  average  energy  per  symbol 

E,  =  a^A^T.  (4.37) 
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The  average  power  per  symbol  P  is 


P  = 


T 


=  a^A^ 


(4.38) 


The  signal  power  at  the  input  to  either  of  the  two  AWGN  ehannel  bloeks  is  the  sum 
power  of  two  multipath  channel  blocks.  Each  of  the  channel  blocks  has  equal  average 
signal  power.  Then  the  total  average  signal  power  P^  at  either  of  the  two  AWGN  blocks 
is 


P,=2P  =  2(7^  . 


(4.39) 


The  total  average  signal  power  P^  must  equal  1 .0  Watt.  The  amplitude  A  is  equal  to  1 
Volt.  Using  this  information  in  Equation  (4.39),  we  find 


l  =  2o-" 

2  1  (4.40) 

<j  =—. 

2 


Thus,  the  variances  of  all  the  Gaussian  noises  sources  of  the  multipath  fading  channels 
are  set  to  1/2 . 


2.  Performance  Analysis 

With  the  addition  of  the  multipath  fading  channels,  at  each  receive  antenna  the  re¬ 
ceived  signals  are 


r'  =  +h^2X^^  +  n^ 

=(^h2iX^  +h22X^^  +  n^. 


(4.41) 


Once  again,  incorporating  the  definitions  for  and  from  Chapter  III,  Equation  (4.41) 
becomes 


r 


1 


r 


2 


=  (/2ii  [xj , -x* ]  +  /2i2  [.^2  Ai*])  +  \n\ ,  n\  ] 
=  [Kx  [^1 , -U ]  +  Ki  [u  Ai  ])  +  [nl 


(4.42) 
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The  channel  coefficients  do  not  change  over  the  code  sequence  and  the  noise  acts  on  one 
symbol  period.  After  serial-to-parallel  conversion,  the  signals  in  each  space  diversity  re¬ 
ception  and  time  interval  are  given  by 


and 


:  [r/ ,  rl  ]  =  [/?!  jXi  +  /212X2  +  nj ,  /212X1  -  h, ^xl  +  nj  ]  . 

(4.43) 

,  ^2^  ]  =  [^2  Al  +  ^22  A  +  ’  ^22-**  “  ^2  A2  +  ^2  ]  • 

(4.44) 

After  time  diversity  combination  of  the  received  signal  on  antenna  1 ,  zj  is 


zl=r;+(r^) 

=  +  /212X2  +  n\  +  (/2i2^i  - +  n\ )’ 

=  +/2j2  jxj  +{j\2  +^^5 

and  Zj  is 

=  +  /Zj2X2  +  n\  -  \hi2^l  -J\^X*2+n\\ 

=  [t\i  -^2)^1  +(^12  +^n)A  ~n2  ■ 


(4.45) 


(4.46) 


Similarly,  the  time  diversity  combination  of  the  received  signal  on  antenna  2  for  x^  and 
X2 ,  represented  by  zf  and  z^  ,  respectively,  are 


zf  =  (^^21  +  ^22  )  -^1  +  (^22  “  ^21  )  A  + 

Zj  =(/221  ~h*27)^\  +(^22  +^2i)a  +  ~  ■ 


(4.47) 


After  space  diversity  combining,  Zj  is 


1  ,  2 

Zj  =  Zj  +  Zj 

=  (/Zjj  +  /z*2  +  /221  +  ^22  )-^l  +  {Ki  ~  K\  +  ^22  “  ^21  )  -^2  +  K  +  ^2  +  ^ 


(4.48) 
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(4.49) 


and  Zj  is 


Z2  =  z*  +  Z2 


[[hn-h 


4-/2  —  h 

12  ^  '^21  '^22 


)xj+(/2, 


12  ^1  ^22  ^21 


1  1*  2  2* 
+  n\  -  ^2  +  . 


Invoking  the  symmetry  of  reeeiver,  we  obtain  the  overall  probability  of  error 
by  using  Zj .  The  demodulator  effeetively  takes  the  real  part  of  Zj .  Onee  again,  let  z  be 
the  real  part  of  Zj .  Then  z  is 

z  =  Re[zj] 

=  ( Re  [A;  j  ]  +  Re  \j1l2  ]  +  Re  [/z2i  ]  +  Re  A22  ]]  Re  [xj  ] 

/  r  n  r  (4-50; 

+  (Re  [^12  ]  -  Re  [hi,  J  +  Re  ]  -  Re  [hi,  ])  Re  [x2  ] 

+  Re  \ji\  j  +  Re [^2*  ]  +  \jh  ]  +  Re  ]  • 

The  real  parts  of  any  h-j  are  IID,  Gaussian  random  variables  with  zero-mean  and 

equal  varianee.  Their  sum  is  also  a  Gaussian  random  variable  with  mean  of  zero  and 
varianee  equal  to  the  sum  of  the  varianees  [5].  As  already  diseussed,  the  real  parts  of  the 
AWGN  eomponents  are  IID  zero-mean  with  equal  varianee.  Their  sum  is  also  a  Gaus¬ 
sian  random  variable  with  zero-mean  and  varianee  equal  to  the  sum  of  the  varianees.  Let 
the  sum  of  ehannel  eoeffieients  multiplied  with  Xj  be  .  Let  the  sum  of  the  ehannel  eo- 

effieients  multiplied  with  X2  be  \  .  Let  the  sum  of  the  noise  eomponents  be  represented 
by  n  .  The  real  parts  of  Xj  and  X2  are  equal  to  Xj  and  X2 .  Equation  (4.50)  beeomes 


z  =  A^Xj  +  +  n  . 


(4.51) 


The  overall  probability  of  bit  error,  ,  ean  be  determined  by  taking  the  sum  of 
probability  of  error  over  all  eombinations  of  Xj  and  Xj  [5].  If  A'  =  a/ 42  ,  then  the  real 
part  of  Xj  and  X2  over  all  eombinations  of  and  ^2  are 
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(4.52) 


b2 

Xi 

X2 

0 

0 

A' 

A' 

0 

1 

A' 

-A' 

1 

0 

-A' 

A' 

1 

1 

-A' 

-A'. 

Now  is  given  by 

=  Pr  {error  |  -  00}  {^1^2  =  00} 

+  Pr{error  |  bjj2  =  0l}Pr{/)i/)2  =  01} 

+  Pr  { error  I /)jZ)2  =  1 0}  Pr  { =  1 0} 

+  Pr  {error  |  b^b^  =  1 1}  Pr  [bjj^  =11} 

The  probability  b^  and  b^  being  a  bit  zero  or  one  are  equally  likely;  therefore, 

Pr(/)j/)2  =  00)  =  Vx{bp2  =  01)  =  Pr(Z)j/)2  =  10)  =  Pr(Z)j/)2  =  1 1)  =  1/4  and  Equation  (4.53)  re- 
duees  to 


P^  =  —(Pr{ error  |  bjj^  =  00}  +  Pr  {error  |  b^b^  =01} 

4  (4.54) 

+  Pr  {error  |  bp^  =  1 0}  +  Pr  {error  |  b^b^  =  1 1}) 

In  terms  of  and  ,  Equation  (4.54)  beeomes 

=  /(']+ Pr[Z  <  0 1 Z,  =  /(',  jr,  =  -A'] 

+  Pr[Z  >  0  \  X^  =  -A',X2  =  A']  +  Pr[z  >  0  \  X^  =  X2  =  -A']). 


Coneentrating  on  the  first  term  of  Equation  (4.55)  and  incorporating  this  informa¬ 
tion  in  Equation  (4.51),  we  get 


z  =  h^A'  +  hi^A'  +  n  . 


(4.56) 


Equation  (4.56)  represents  a  linear  transformation  of  three  independent  Gaussian 
random  variables  and  n  .  The  linear  transformation  of  Gaussian  random  variables 

remains  Gaussian  [5].  Therefore,  z  is  a  Gaussian  random  variable  with  mean  equal  to 
the  sum  of  the  means  and  variance  equal  to  the  sum  of  the  variances.  Each  of  and 
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n  have  zero-means.  Thus,  z  has  a  zero-mean.  The  first  term  of  Equation  (4.55)  is  given 
by  [5] 

Pr[Z<0,Xi  =^2  =^']  =  Pr[/2^^'  +  /2,^'  +  n<0]  =  |.  (4.57) 


Each  of  the  three  remaining  terms  in  Equation  (4.55)  gives  the  same  result;  thus,  the 
overall  probability  of  bit  error  is 


Pe- 


1 

4 


1111 

- 1 - 1 - 1 - 

2  2  2  2 


1 

2  ■ 


(4.58) 


Therefore,  the  probability  of  bit  error,  ,  for  the  system  is  0.5,  regardless  of  the  of  the 
signal  or  the  power  of  the  noise.  This  was  confirmed  by  simulation  for  increasing  aver¬ 
age  symbol  energy  to  noise  ratio  EjN^  . 


3.  Performance  Simulation  Results 

The  simulation  was  conducted  with  increasing  EjN^  .  The  period  T  was  set  to 

0.001  seconds.  The  simulation  run  time  was  set  at  100.0  seconds.  Thus,  the  total  num¬ 
ber  of  bits  in  a  single  simulation  run  were  100,000. 

The  simulated  probability  of  error  for  the  MIMO  system  is  plotted  in  Eigure  17. 
The  results  are  plotted  with  the  theoretical  results  obtained  in  Equation  (4.58)  In  addi¬ 
tion,  the  theoretical  probability  of  error  for  a  baseband  equivalent  SISO  BPSK  system  in 
a  multipath  channel  derived  in  Chapter  II  is  plotted  for  comparison.  This  is  given  by 

Pe  = 


1- 


r, 

\+f, 


(4.59) 


where  is  given  by 


r,= 


Eqcj^ 

N. 


(4.60) 


Recall  from  Equation  (4.40)  that  cr^  =  1/2  thus.  Equation  (4.60)  becomes 
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(4.61) 


The  simulated  results  follow  the  theoretical  results  very  well.  However,  they  do  not  give 
the  performance  expected  by  a  baseband  equivalent  SISO  BPSK  in  a  multipath  channel. 


Figure  17.  Results  in  Multipath  Fading  Channel. 

The  conclusion,  of  course,  is  that  the  receiver  as  designed  in  Chapter  111  and  simu¬ 
lated  in  Figure  12  cannot  communicate  in  a  multipath  fading  channel.  Thus,  improve¬ 
ment  to  the  receiver  design  must  be  incorporated  in  order  for  the  system  to  function  in  a 
multipath  channel.  In  the  next  section,  the  required  improvement  is  incorporated  in  the 
design;  the  MRC  receiver  with  CSl. 
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D,  SIMULATION  IN  A  MULTIPATH  CHANNEL  WITH  MRC  DESIGN 

An  MRC  receiver  incorporates  CSI  [10].  The  incorporation  of  CSI  follows  that 
found  in  Reference  [13].  Nevertheless,  before  discussing  the  design  of  the  receiver,  the 
received  signals  are  examined  more  closely.  Investigation  of  Equations  (4.43)  and  (4.44) 
once  again  gives 


r/  =  +/212X2  +n\ 

=  /Zj2Xj  -/ZjjX2  +^2 

+  ^22X2  +  nf 


2  r  *  j  *  ,  2 

=  h22Xi  -/22A2  +^2 


(4.62) 


Next,  we  take  the  complex  conjugate  of  and  ,  and  Equation  (4.62)  becomes 


q'  +^12-^2  +^1 

(^2* )  ~  ~  Ki^2 + 

+  A22X2  +  nf 

(2 \*  7  *  7  *  2* 

^2  )  “  ^22'^1  “  ^21-^2  ^2  • 


(4.63) 


We  would  be  well  served  to  view  Equation  (4.63)  as  describing  a  2-input,  4-output 
MIMO  system.  This  equation  set  can  be  represented  in  matrix  form  as 


where  r  is 


X  is 


r  =  Hx  -I-  n , 


X  = 


X, 

X2 


(4.64) 


(4.65) 


(4.66) 
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H  is 


\\  \i 


“11 

“12 

^21 

h 

“22 

hn 

-hn 

hn 

-hn 

and  n  is 


n  = 


(4.67) 


(4.68) 


Note  that  the  second  and  third  lines  in  Equation  (4.63)  were  swapped  in  the  definitions 
described  in  Equations  (4.66),  (4.67)  and  (4.68). 

Erom  Chapter  II  and  Reference  [I],  the  capacity  is  dependent  on  matrix  Q  given 
in  Equation  (2.74)  by 


Q  =  H’^H . 


Using  the  definition  of  H  in  Equation  (4.67),  we  get 


Q  =  H*  H 

_  \\  ^21  \l 

which  reduces  to 


^22 

-h 


21 


/ill  ^2 

/It,  h 


'22 

-k 


11 


a;  -h 


21 


(4.69) 


(4.70) 


(4.71) 


As  discussed  in  Chapter  II  and  Reference  [I],  this  is  a  form  of  matrix  Q  that 
gives  a  linear  increase  in  capacity  with  size  of  Q  ,  in  this  case  two.  If  CSI  is  used 
appropriately,  then  this  capacity  can  be  approached.  In  the  receiver,  if  channel 
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priately,  then  this  capacity  can  be  approached.  In  the  receiver,  if  channel  information  is 
used  according  to  ,  then  this  can  be  achieved. 

It  seems  extraordinary  that  the  matrix  H  has  this  property.  However,  it  was  the 
code  and  the  reception  of  the  code  that  forced  this  property  [13].  Referring  to  Chapter 
III,  we  see  that  the  code  matrix  X  had  this  same  property.  Thus,  in  receiving  this  code, 
the  channel  matrix  H  has  this  property  regardless  of  h.j ,  the  channel  realization  over  the 

code  interval  [13]. 

The  use  of  CSI  must  be  incorporated  according  to  .  Given  the  received  ma¬ 
trix  r ,  let  a  signal  matrix  z  be  defined  by 

z  =  H’^r  .  (4.72) 

This  is  consistent  with  Reference  [13].  Thus,  the  implementation  of  CSI  is  accomplished 
according  to  Equation  (4.72). 


The  expansion  of  Equation  (4.72)  gives 


^11  ^21  ^12  ^22 

'X 

hi,  a;  -h,,  -h,,_ 

0 

- 1 

* 

+  ki  (c )  +  Ki  (c 

Ki^x  ^Ki^x  -Kx[4) 

^21  (C 

(4.73) 


1.  Simulation  Parameters 

Equation  (4.73)  represents  the  rules  for  incorporating  the  CSI.  The  modified  re¬ 
ceiver  that  incorporates  the  use  of  CSI  is  described  in  Eigure  18.  The  modified  receiver 
is  a  MRC  receiver.  The  simulation  of  the  MRC  receiver  is  described  in  Eigure  19.  The 
channel  parameters  ^22  taken  directly  from  the  appropriate  multipath 

fading  block.  This  provides  perfect  channel  knowledge. 
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Figure  18.  MIMO  Receiver  with  CSI  implementation. 


Figure  19.  Multipath  Simulation  with  MIMO  MRC  receiver. 


In  order  to  coincide  with  the  correct  symbols  at  the  multipliers,  the  channel  coefficients 
are  delayed  by  2T  seconds  since  the  buffers  in  the  receiver  induces  this  delay.  The  re¬ 
mainder  of  the  simulation  parameters  remain  the  same  as  in  previous  simulations. 


2,  Performance  Analysis 

The  received  signal  matrix  z  represents  the  signals  Zj  and  Zj  given  by 


z  = 


H*"r, 


(4.74) 


which  expands  to 


z  =  H’^Hx  +  H*^n 


z=Qx  +  H’^n. 

After  matrix  multiplication,  Equation  (4.75)  expands  to 


z  = 


+|^2l|^  +1^12^  +|^22r)-^l  +^2^r  +h*2l^l  +  ^22^f 

||^ll|  "^1^211  '^1^22!  j -^2  ~  ^1^2  ^12^1  ~  ^21^2  ^22^1 


(4.75) 


(4.76) 


Applying  the  symmetry  of  the  receiver,  we  only  have  to  investigate  Zj .  The  result 
of  Equation  (4.76)  equates  to  an  MRC  receiver  with  L  =  4  diversity  receptions  [10].  The 
probability  of  error  of  Zj  was  given  in  Chapter  II,  that  is. 


R  = 


(l-w) 

^L-l 

^L-\  +  t 

V  / 

2 

SI 

/=0 

1  1  ) 

(4.77) 


where  u  is 


u  = 


(4.78) 
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and  y^  is 


Yi  = 


EIct^ 


(4.79) 


The  parameter  E  in  Equation  (4.79)  represents  the  energy  per  symbol  E^ .  The 
energy  per  symbol  is  given  by 


Since  cr^  =  1/2  ,  Equation  (4.79)  becomes 


(4.80) 


(4.81) 


3,  Performance  Simulation  Results 

The  simulation  was  conducted  with  increasing  EjN^  .  As  was  the  case  in  the 

AWGN  simulation,  the  period,  T ,  was  set  to  0.0001  seconds.  The  simulation  run  time 
was  set  at  500.0  seconds  giving  the  total  number  bits  in  a  single  simulation  run  of  five 
million. 

The  simulated  probability  of  error  for  the  MIMO  systems  is  plotted  in  Eigure  20. 
The  results  are  plotted  with  the  theoretical  results  given  by  Equations  (4.77),  (4.78)  and 

(4.81).  As  well,  the  theoretical  probability  of  error  for  a  SISO  BPSK  system  is  plotted 
for  comparison.  Note  that  the  SISO  BPSK  system  has  perfect  knowledge  of  the  CSI. 

The  MIMO  system’s  simulated  results,  with  MRC  receiver,  follow  the  theoretical  results 
very  well.  The  simulated  results  deviate  starting  at  approximately  13.0  dB.  This  is  a  re¬ 
sult  of  the  number  of  bits  per  simulation  run.  Overall,  the  simulation  follows  the  theo¬ 
retical  analysis  of  bit  error  performance.  Also  indicated  in  the  figure  is  the  improvement 
in  bit  error  performance  over  the  single  channel  BPSK  system.  Eor  equal  E^^jN^ ,  there 

is  an  excellent  improvement  in  performance  that  increases  as  Ej^jN^  increases.  This  im¬ 
provement  follows  the  expected  improvement  as  illustrated  in  Chapter  II.  Easily,  the  sys¬ 
tem  achieved  E  =  4  diversity  order;  specifically,  the  sum  of  two  time  diversity  and  two 
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space  diversity  orders.  This  is  consistent  with  the  potential  of  the  Alamouti  space-time 
code  [1,13], 


Figure  20.  Results  MRC  Multipath  fading. 

Found  in  Reference  [1]  are  the  results  for  probability  of  error  for  a  MIMO  system 
using  the  Alamouti  scheme.  They  are  given  in  Figure  21 .  The  published  results  follow 
the  same  assumptions,  i.e.,  the  channels  are  independent  and  perfect  CSI  is  known  [1].  In 
the  figure,  the  graph  with  the  circle  markings  is  the  one  of  primary  interest,  which  corre¬ 
sponds  to  the  simulation  in  this  thesis.  The  graph  indicates  very  similar  results,  particu¬ 
larly  at  =  10  dB ,  where  is  approximately  10 . 


E.  SUMMARY  OF  RESULTS 

In  this  chapter,  the  MIMO  system  designed  in  Chapter  III  was  examined  with 
A  WON  only  as  well  as  the  effects  of  multipath.  The  MIMO  system  performed  as  ex- 


80 


pected  with  only  AWGN  and  achieved  better  performance  than  a  SISO  BPSK  system. 
The  system  was  than  simulated  for  a  multipath  ehannel.  As  designed,  it  was  not  able  to 
eommunicate  through  this  channel.  The  receiver  for  the  original  system  was  converted  to 
an  MRC  receiver  with  CSI.  As  a  result,  the  new  design  was  able  to  eommunieate  suc- 
eessfully  in  the  multipath  ehannel.  The  performanee  of  the  new  design  matehed  expeeted 
theoretieal  performance.  Additionally,  it  also  outperformed  a  SISO  BPSK  system. 


Figure  21.  Published  results  (From  Reference  [1].). 


Furthermore,  the  system’s  performance  compared  well  with  published  results.  The 
MIMO  system  was  able  to  aehieve  T  =  4  diversity  order  with  the  Alamouti  space-time 
code,  whieh  is  eonsistent  with  published  results. 

In  the  next  ehapter,  the  results  of  this  thesis  are  summarized  and  areas  for  follow- 
on  work  are  presented. 
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V.  CONCLUSION 


The  goal  of  this  thesis  was  to  develop  a  MIMO  system  with  space-time  coding 
with  an  eye  towards  military  application.  A  MIMO  system  was  designed  using  the 
Alamouti  space-time  code.  It  was  simulated  in  Matlab  with  Simulink  and  was  tested  in  a 
progressive  manner,  first  in  an  AWGN  channel  and  then  in  a  multipath  fading  channel 
with  AWGN.  The  performance  of  the  system  was  analyzed  and  compared  with  simulated 
results,  as  well  as  published  results  found  in  Reference  [1]. 

A.  RESULTS 

The  designed  MIMO  system  performed  in  a  manner  consistent  with  theoretical 
analysis  in  the  AWGN  channel.  The  simulated  results  compared  very  well  with  the  theo¬ 
retical  performance  analysis.  However,  the  original  design  failed  in  a  multipath  channel 
with  AWGN.  The  design  did  not  incorporate  CSI  and  a  MRC  receiver.  Nevertheless, 
performance  analysis  was  derived  and  simulated  results  reflected  the  anticipated  theoreti¬ 
cal  results.  With  this  failure,  the  original  receiver  design  was  changed  to  an  MRC  re¬ 
ceiver  incorporating  CSI.  Consequently,  the  new  design  was  able  to  communicate  suc¬ 
cessfully  in  the  multipath  channel.  The  performance  of  the  new  design  matched  expected 
theoretical  performance,  as  well  as  published  results.  Finally,  the  designed  MIMO  sys¬ 
tem  with  the  MRC  receiver  was  able  to  achieve  full  diversity  order  with  the  Alamouti 
space-time  code  which  was  consistent  with  published  results.  There  is  improved  bit  error 
performance  in  the  MIMO  system  over  the  SISO  system  and  capacity  gain  which  can  be 
exploitable  for  military  application. 


B,  FOLLOW-ON  WORK 

There  are  three  areas  identified  for  follow-on  work.  First,  the  designed  MIMO 
receiver  can  be  converted  into  a  maximum-likelihood  decoding  receiver.  This  design 
holds  promise  of  simplicity  over  the  original  design.  The  performance  analysis  can  be 
investigated  and  simulated. 
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Next,  higher  order  spaee-time  eodes  ean  be  investigated.  Although  the  Alamouti 
space-time  code  can  achieve  full  diversity  gain,  it  is  unable  to  achieve  coding  gain  as 
discussed  in  Chapter  II  without  bandwidth  expansion.  Higher  order  codes  can  be  exam¬ 
ined  and  tested  by  augmenting  the  original  design  or  using  a  maximum-likelihood  decod¬ 
ing  receiver.  The  incorporation  of  ECC  can  be  investigated  in  conjunction  with  the 
modulation  scheme  in  order  to  achieve  coding  gain. 

Lastly,  the  MIMO  system  is  already  using  CSI  at  the  receiver.  If  this  information 
about  the  channel  can  be  extracted  at  the  receiver,  then  it  can  be  sent  to  the  transmitter. 
The  transmitter  can  use  CSI  to  more  efficiently  transmit  power  to  the  receiver  through 
channels  that  are  less  faded  than  others.  Techniques  that  control  the  transmitter  can  be 
investigated  and  performance  of  the  system  analyzed. 
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APPENDIX  A.  LIST  OF  ACRONYMS 


AWGN 

Additive  White  Gaussian  Noise 

BPSK 

Binary  Phase-Shift  Keying 

CSI 

Channel  State  Information 

ECC 

Error  Control  Coding 

IID 

Independent  Identieally  Distributed 

MIMO 

Multiple-Input-Multiple-Output 

MRC 

Maximal-ratio  Combining 

SISO 

Single-Input-Single-Output 

SNR 

Signal-to-Noise  Ratio 

SVD 

Singular  Value  Deeomposition 
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APPENDIX  B.  MATLAB  SIMULINK  BLOCK  PARAMETERS 


The  block  parameters  shown  in  this  appendix  give  the  parameters  and  a  brief  de¬ 
scription  of  the  Matlab  blocks  used  in  the  simulations  in  this  thesis. 
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Block  Parameters:  Interlacerl 


-Interlacer  (mask)  (link) - 

Combine  the  elements  of  the  input  signals  to  generate  the  output  signal. 
The  elements  of  the  first  input  signal  become  the  odd-numbered  elements 
of  the  output  signal,  while  the  elements  of  the  second  input  signal  become 
the  even-numbered  elements  of  the  output  signal. 

The  inputs  can  be  either  scalars  or  frame-based  column  vectors. 


OK 

Cancel 

Help  1 

Apply  j 
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Block  Parameters:  AWGN  Channell 


x] 

-AWGN  Channel  (mask)  (link) - 

Add  white  Gaussian  noise  to  the  input  signal.  The  input  and  output  signals 
can  be  real  or  complex.  This  block  supports  multichannel  input  and  output 
signals  as  well  as  frame-based  processing. 


When  using  either  of  the  variance  modes  with  complex  inputs,  the 
variance  values  are  equally  divided  among  the  real  and  imaginary 
components  of  the  input  signal. 


OK  I  Cancel  |  Help  |  Apply  | 
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